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ABSTRACT 


The  problem  of  efficiently  packing  many  mobile  user  signals  into  the  avail- 
able UHF  uplink  bandwidth  of  an  advanced  communications  satellite  is  considered. 
It  is  assumed  that  the  user  terminals  are  mutually  unsynchronized  in  timing, 
frequency  and  transmitter  power  control,  and  that  the  satellite  performs  the 
processing  functions  of  tracking  the  symbol  timing  and  Doppler  shift  of  each 
user,  demodulating  the  uplink  signals,  and  reorganizing  the  information  for 
the  downlink.  The  basic  approach  is  FDMA  with  QPSK  modulated  data  streams 
individually  modified  by  spectral  shaping  using  data  windows.  BPSK  and  MSK 

modulation  formats  were  also  investigated. 

From  extensive  computer  simulations  it  is  concluded  that  with  relatively 

little  degradation  20  unsynchronized  16  kbps  satellite  users  can  be  packed 

into  a 500  kHz  uplink  bandwidth  that  is  shared  with  terrestrial  LOS  allocations. 

Alternatively,  7 unsynchronized  2400  bps  users  can  occupy  a 25  kHz  bandwidth  or 

133  such  users  can  be  packed  into  a 500  kHz  bandwidth.  Signal  strengths  may 

vary  by  as  much  as  10  to  15  dB  at  the  satellite  without  the  need  for  transmitter 

power  control.  Peak-power-limited  QPSK  with  a rate  1/2  convolutional  code  is 

_3 

slightly  inferior  to  uncoded  MSK  for  bit  error  rates  larger  than  10  . Spec- 

trally shaped  and  coded  QPSK,  however,  is  needed  for  smaller  error  rates  and 
will  be  more  attractive  if  the  modulator  is  not  peak-power-limited. 
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I.  SUMMARY 

The  subject  of  this  report  is  modulation  techniques  for  bandwidth  conser- 
vation utilizing  an  advanced  processing  satellite  envisioned  for  the  1980's. 

1.1  Introduction 

There  is  a need  to  provide  efficient  service  to  many  mobile  Ultra  High 
Frequency  (UHF)  terminals  operating  in  a frequency  division  multiple  access  (FDMA) 

uplink  mode  [1].  These  terminals  may  be  located  on  submarines,  ships,  air- 
craft, etc.,  operating  anywhere  on  the  surface  of  the  earth.  The  terminals 
will  operate  at  a variety  of  communication  data  rates  [2].  Three  representa- 
tive data  rates  are  75  bits  per  second  (bps)  (teletype),  2400  bps  (data  or 
vocoded  speech)  and  16  kbps  (continuously  variable  slope  delta  (CVSD)  modula- 
tion voice  [3,  4],  which  has  a particularly  simple  implementation).  Since  the 
terminals  will  vary  considerably  in  effective  radiated  power  (ERP)  the  received 
signal  powers  at  the  satellite  can  differ  by  as  much  as  15  or  20  dB. 

Because  a large  number  of  terminals  is  anticipated,  the  overall  system  cost 
dictates  that  the  user  terminals  be  relatively  simple.  Transmitter  power  con- 
trol, symbol  timing  and  Doppler  correction  all  contribute  to  terminal  complexity. 
It  is  generally  desirable  to  shift  the  burden  for  realizing  these  functions  to 
the  satellite  as  much  as  possible  in  order  to  hold  down  the  cost  of  a typical 
terminal. 

Another  concern  is  the  problem  of  frequency  allocation  in  the  government 
UHF  band  (225  to  400  MHz).  There  are  thousands  of  line-of-sight  (LOS)  users 
allocated  at  100  kHz  centers.  In  anticipation  of  more  closely  spaced  alloca- 
tions, there  are  a large  number  of  terminals  capable  of  tuning  to  50  kHz  centers 
and  a smaller  but  growing  number  tunable  to  25  kHz  centers  [5-7].  It  is 
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unlikely  that  more  than  a few  hundred  kHz  of  contiguous  bandwidth  would  be 

available  for  the  allocation  of  satellite  users.  Thus,  the  uplink  will  probably 

be  composed  of  a collection  of  sub-bands,  each  one  of  which  might  contain  a 

* 

number  of  satellite  users  and  perhaps  some  LOS  users  as  well. 

If  a degree  of  anti-jam  (AJ)  protection  is  desired,  these  sub-bands  could 
also  be  frequency  hopped  at  a modest  rate,  e.g.,  75  hops  per  sec  to  be  com- 
patible with  the  teletype  rate.  The  problem  of  reacquiring  phase  at  the  re- 
ceiver with  each  frequency  hop,  a difficulty  that  is  present  with  coherent 
baseband  modulation  schemes,  is  not  addressed  in  this  report. 

An  FDMA  uplink  mode  is  an  attractive  candidate  under  these  circumstances. 
Although  certain  frequency  bands  could  actually  contain  time  division  multiple 
access  (TDMA)  traffic,  this  possibility  will  not  be  developed  here.  Two  pro- 
blems with  TDMA  are  that  the  weaker  terminals  may  not  be  able  to  support  the 
burst  rates  required,  and  the  terminals  would  have  to  be  synchronized. 

Looking  toward  the  future  in  satellite  communications  it  is  reasonable 
to  expect  significant  progress  in  on-board  processing  capabilities  that  would 
be  reflected  in  upgraded  service  in  several  areas.  Ten  years  from  now  a gen- 
eral purpose  UHF  satellite  could  provide  at  least  500  kbps  of  throughput  data 
rate,  which  is  an  order  of  magnitude  more  capacity  than  will  be  available  in 
the  near  term.  Advanced  modulation  techniques  would  permit  the  more  efficient 
use  of  the  crowded  UHF  band.  Separation  of  the  uplinks  and  downlinks  would 
result  in  greater  system  versatility  emcompassing  such  features  as  more  cross- 

* 

Naturally,  it  is  important  that  the  satellite  users  and  the  LOS  users  do 
not  interfere  with  one  another.  This  issue  will  be  addressed  by  the  spectral 
shaping  of  the  satellite  user  signals. 
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band  operation  between  UHF  and  other  frequency  bands,  selective  treatment  of 
different  coverage  areas  on  the  earth,  and  uplink  AJ  protection  at  UHF. 

Thus,  the  satellite  itself  is  presumed  to  be  an  advanced  processing 
satellite  capable  of  demodulating  the  uplink  signals,  reorganizing  the  demodu- 
lated data  for  the  downlink  and  remodulating  the  information  according  to  the 
destinations  specified.  Such  a satellite  might  employ  multiple  beam  antennas, 
perform  AJ  antenna  array  nulling  and  frequency  dehopping,  participate  in  the 
tracking  and  correction  of  symbol  timing  and  Doppler  shifts  for  uplink  users, 
and  execute  the  multifarious  control  functions  necessary  for  such  versatile 
operation. 

For  the  purpose  of  this  report  the  satellite  downlink  may  be  thought  of 
as  either  frequency  division  multiplexed  (FDM)  or  time  division  multiplexed  (TDM) . 
The  reasons  for  preferring  one  or  the  other  downlink  mode  are  beyond  the  scope 
of  this  report  and  will  not  be  discussed. 

The  fundamental  objective  of  this  report  is  to  establish  a means  for  ac- 
complishing the  close  frequency  packing  of  more  than  the  usual  number  of  satel- 
lite users  in  a contiguous  section  of  uplink  bandwidth. 

With  the  hope  of  achieving  a minimum  of  terminal  and  satellite  processing, 
this  close  packing  could  be  attempted  for  unsynchronized  users  (for  terminal 
simplicity)  by  carefully  selecting  an  appropriate  modulation  scheme.  The  de- 
gree of  packing  would  be  limited  by  the  center  frequency  separations  required 
to  reduce  interuser  interference  to  a tolerable  level. 

If  the  users  were  synchronized  the  modulated  waveforms  arriving  at  the 
satellite  could  be  made  orthogonal  and  the  users  could  indeed  be  very  closely 
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packed  in  frequency.  Furthermore,  there  would  be  no  power  control  problem 
because  of  the  signal  orthogonality.  On  the  other  hand,  synchronization  in- 
creases terminal  complexity  and  systems  dependent  on  sync  may  be  fragile. 

Hence,  it  is  desirable  to  consider  the  efficacy  of  the  unsynchronized  case. 

With  unsynchronized  users  there  will  generally  be  some  crosstalk  (CT) 
among  the  various  received  signals  in  the  satellite  demodulation  channels.  As 
a first  step  in  attaining  the  close  packing  objective,  the  worst  case  CT  be- 
tween just  two  uncoordinated  user  signals  has  been  calculated  for  eight  differ- 
ent types  of  simple  modems,  including  those  utilizing  binary  phase  shift 
keying  (BPSK) , a few  variants  of  quadriphase  shift  keying  (QPSK) , and  continu- 
ous phase  frequency  shift  keying  (CPFSK) . These  results  are  interesting  in 
themselves  but  because  of  the  rather  restricted  scope  of  the  investigation  it 
was  decided  to  treat  that  subject  in  a separate  but  companion  report  [8].  How- 
ever, the  major  conclusion  of  that  study  will  be  included  in  the  next  section 
of  this  report. 

In  this  report  the  principal  concern  is  with  the  more  general  situation 
where  additional  terminal  and/or  satellite  processing  is  permitted.  In  parti- 
cular, user  signals  may  be  either  synchronized  or  unsynchronized  in  time 
and/or  frequency.  Spectral  shaping  of  the  transmitted  and/or  received  signals 
will  allow  users  to  be  closely  spaced.  The  shaping  is  accomplished  by  digital 
filtering  at  either  the  terminal  or  the  satellite  or  both. 

The  general  system  block  diagram  indicating  the  approach  taken  in  this 
report  is  given  in  Fig.  1.1.  Since  most  of  the  processing  for  bandwidth  con- 
servation is  best  accomplished  digitally,  and  because  the  mathematics  of  analog 
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and  digital  processing  are  basically  equivalent  [9],  very  little  will  be  said 

about  the  analog  or  digital  to  analog  (D/A)  portions  of  the  system.  The  analog 

to  digital  (A/D)  portion  will  be  discussed  briefly  in  Section  5.1.  Specific 

options  for  this  mix  between  analog  and  digital  satellite  processing,  multiple 

front  ends,  intermediate  frequencies  (IF!s),  etc.,  can  have  a considerable 
impact  on  satellite  design  but  these  issues  will  be  left  for  future  study. 

As  seen  in  Fig.  1.1  just  one  contiguous  section  of  user  bandwidth  is  as- 
sumed. At  some  point  in  the  satellite  a digital  baseband  signal  is  obtained 
and  fed  to  a number  of  replicated  user  demodulation  channels,  one  channel  for 
each  user  to  be  demodulated.  Although  the  primary  focus  of  this  report  is  on 
a typical  demodulation  channel  there  will  be  an  attempt  to  deal  with  some  of 
the  larger  satellite  processing  questions  involving  synchronized  versus  unsyn- 
chronized users,  timing  and  Doppler  tracking,  mixed  data  rates,  etc. 

In  this  report,  it  is  assumed  that  the  satellite  tracks  the  symbol  timing 
and  Doppler  of  each  user  it  demodulates  without  elaborating  on  how  this  might 
be  done. 

The  essence  of  the  windowing  approach  of  this  report  will  now  be  described. 
As  depicted  in  Fig.  1.1  the  filtering  for  spectral  shaping  is  applied  in  the 
time  domain  by  amplitude  weighting  (windowing)  the  signals  at  the  transmitter 
and  satellite  with  data  windows.  The  desired  weighting  coefficients  can  be 
generated  easily  with  digital  read  only  memories  (ROM’s)  and  D/A  converters. 

For  the  problem  considered  here,  this  is  more  efficient  than  convolutional 
filtering  [10-12].  In  addition,  data  windows  are  finite  in  extent  which 
enables  the  windows  to  be  applied  without  introducing  intersymbol  interfer- 
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ence  (ISI) . ISI,  when  introduced  by  windowing  over  more  than  one  information 
symbol  as  shown  in  Fig.  1.1,  can  be  easily  mitigated.  As  will  be  seen  the  ISI 
is  specified  by  the  crosscorrelation  function  of  the  satellite  and  transmitter 
windows . 

Windowing  at  the  terminal  controls  the  spectrum  of  the  transmitted  signal 
so  that  there  is  less  spectral  splatter  from  one  satellite  user  to  another. 
Transmitter  windowing  will  also  reduce  the  crosstalk  (CT)  between  a satellite 
user  and  a line  of  sight  (LOS)  user  in  the  vicinity  of  the  transmitting  satel- 
lite user.  Windowing  at  the  satellite  reduces  out  of  band  interference  and 
restricts  the  demodulation  bandwidth  for  a given  satellite  user.  If  the  spec- 
tral window  at  the  satellite  is  sufficiently  narrow,  it  will  also  reject  the 
signal  from  an  LOS  user  at  the  band  edge.  In  general,  the  longer  the  data 
windows  the  more  effective  the  spectral  shaping  will  be  and  the  less  the  CT. 

On  the  other  hand,  longer  windows  generally  mean  increased  ISI.  Transmitter 
windowing  may  also  imply  non-constant  envelope  signaling  which  for  peak-power- 
limited  transmitters  means  a received  signal  energy  penalty  with  respect  to 
channel  noise.  So  in  this  model  of  the  FDMA  environment,  there  are  four  sources 
of  error  — the  usual  channel  noise,  assumed  to  be  additive  white  Gaussian 
noise  (AWGN) , CT  among  users,  ISI,  and  LOS  interference.  Clearly,  there  will 
be  system  design  tradeoffs  among  these  four  sources  of  interference. 

There  are  many  ways  to  mitigate  ISI  in  the  presence  of  AWGN.  In  relation 
to  most  situations  covered  in  the  literature  [13  - 23],  the  ISI  problem  at  hand 
is  rather  a mild  one  since  the  ISI  is  deliberately  introduced  (hence,  known). 
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of  limited  extent  and  t ime- invar ian t . As  will  be  shown  in  this  report,  the  ISI 

that  arises  in  controlling  CT  from  closely  packed  users  can  be  effectively 

eliminated  with  a simple  digital  filter.  As  indicated  in  Fig.  1.1  the  filter 

can  be  designed  to  minimize  the  mean  square  error  (MSE)  between  the  input  to 

the  decision  device  in  the  satellite  and  the  actual  symbol  transmitted.  The 

resulting  filter  is  suboptimum  in  the  sense  that  the  probability  of  a symbol 

error  is  not  actually  minimized.  However,  it  is  also  shown  through  computer 

-3 

simulation  that  the  filter  performs  nearly  identically,  at  least  down  to  10 
bit  error  rates,  to  the  more  complex  Viterbi  algorithm  implementation  which 

realizes  a maximum  likelihood  detection  on  the  entire  transmitted  sequence. 

Being  identical  in  design  the  ISI  filters  and  decision  devices  could  be 
either  replicated  or  time  shared  among  user  demodulation  channels. 

Presumably  some  combination  of  partially  coherent  [24  -28]  or  differentially 
coherent  [29,  30]  detection  will  be  employed  in  an  actual  implementation.  The 
theoretical  and  experimental  degradations  from  coherent  performance  in  AWGN 
are  well  established  [31,  32]  and  can  be  used  in  a more  refined  analysis.  For 
simplicity  and  concreteness  only  coherent  matched  filter  hard  decision  re- 
ceivers [33]  are  considered  in  this  report. 

The  baseline  modulation  scheme  employed  is  QPSK  [34],  although  some  of 
the  computer  simulated  results  obtained  may  apply  equally  well  to  offset  QPSK 

[35,  36]  and  M-ary  PSK.  Simulation  has  also  been  performed  for  BPSK  and  mini- 

* 

mum  shift  keying  (MSK)  [37  - 39],  Some  related  work  dealing  with  up  to  two 
interfering  QPSK  user  signals  is  found  in  [40,  41] . 

* 

Binary  continuous  phase  FSK  with  a frequency  deviation  ratio  of  1/2. 
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Coding  issues  are  secondary  to  modulation  techniques  for  bandwidth  con- 
servation and  are  not  considered  with  much  detail  in  this  report.  However,  it 
is  recognized  that  coding  may  be  necessary  to  achieve  the  required  system  per- 
formance especially  when  considering  all  types  of  interference.  Convolutional 
coding  [42]  can  be  a particularly  efficient  means  of  improving  performance 
with  a reasonable  increase  in  system  complexity.  At  present,  coding  is  usually 
considered  to  be  end-to-end  in  nature,  i.e.,  no  uplink  decoding  or  re-encoding 
for  the  downlink  in  the  satellite.  For  convenience  in  evaluating  the  uplink 
performance  with  coding,  decoding  is  assumed  to  take  place  in  the  satellite, 

or  alternatively,  downlink  noise  is  neglected.  This  analytical  device  affects 
the  performance  level  only  slightly  and  has  very  little  bearing  on  the  conclu- 
sions of  this  report  which  are  presented  in  the  next  section. 

The  rest  of  the  report  is  organized  as  follows.  Section  II  contains  a 
description  of  the  mathematical  model  of  the  system  and  an  analysis  of  ways  to 
mitigate  intersyrabol  interference  introduced  by  the  spectral  shaping.  Section  III 
discusses  data  windows  which  produce  the  spectral  shaping  for  controlling  cross- 
talk among  closely  packed  user  signals.  Simulated  results  of  various  system 
designs  are  presented  in  Section  IV.  Finally,  Section  V discusses  some  imple- 
mentation issues  bearing  on  system  complexity.  A glossary  and  a representative 

list  of  references  are  provided  in  this  report  to  assist  the  interested 
reader. 

1.2  Conclusions 

The  overall  results  of  this  report  and  its  companion  [8 ] are  discussed  in 
this  section.  The  major  conclusions  are  capitalized  for  emphasis. 
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The  principal  conclusion  of  the  companion  report  where  crosstalk  (CT)  among 

users  is  combatted  simply  by  sufficient  center  frequency  separation  of  the 

modulated  waveforms  with  no  in-band  filtering  (See  Table  1.1)  is  that: 

WITH  NO  SPECIAL  PROCESSING,  I.E.,  NO  WINDOWING, 

MANY  MORE  CONTINUOUS  PHASE  FSK  USERS  CAN  OCCUPY 
A GIVEN  BANDWIDTH  THAN  EITHER  BPSK  OR  QPSK  USERS. 

A worst  case  analysis  assuming  hard  decision  coherent  receivers  suggests  that 
three  to  five  times  as  many  minimum  shift  keyed  (MSK)  waveforms  as  QPSK  wave- 
forms within  a specified  bandwidth  are  possible  at  a tolerable  level  of  CT 
interference.  The  bandwidth  advantage  of  MSK  over  the  variants  of  QPSK  was 
shown  to  increase  with  the  power  level  of  the  interfering  waveforms. 

At  this  point  a brief  explanation  of  this  behavior  is  worthwhile.  An  MSK 
signal  is  a continuous  waveform  that  can  be  detected  with  the  same  bit  error 

probability  in  AWGN  as  BPSK  or  QPSK.  The  phase  continuity  of  MSK  results  in  a 
very  compact  signal  spectrum  with  a rapidly  decreasing  sidelobe  level  as  a func- 
tion of  frequency.  Although  the  bandwidth  between  the  principal  spectral  nulls 
of  MSK  is  50%  larger  than  that  of  QPSK  for  random  data  at  the  same  rate,  the 
fraction  of  out-of-band  energy  is  generally  much  smaller  for  MSK  than  for  QPSK. 
It  has  been  determined  that  the  crosstalk  decreases  only  linearly  with  center 
frequency  separation  with  M-ary  PSK  while  with  continuous  phase  FSK  the  CT 
dependence  varies  as  the  inverse  square  of  the  frequency  separation.  This  is 
consistent  with  the  asymptotic  behavior  of  the  various  spectra.  The  reader  is 
referred  to  the  companion  report  [8]  for  more  details. 

Later  in  this  section  the  number  of  user  signals  possible  in  a given  band- 
width, and  at  a certain  data  rate  and  performance  level,  both  with  a without 
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windowing,  will  be  estimated.  Naturally,  a greater  density  of  waveforms  will 
be  possible  using  the  special  windowing  techniques  of  this  report.  Subsequent 
discussion  will  focus  on  the  conclusions  that  arose  from  this  additional  sig- 
nal processing  for  bandwidth  conservation. 

In  this  discussion  the  modulated  waveforms  transmitted  by  the  satellite 
users  are  nominally  identical  in  form  and  spectrally  located  at  closely  packed 
but  distinct  center  frequencies.  Either  QPSK  or  MSK  modulation  is  assumed  with 
a user  data  rate  of  16  kbps  and  a center  frequency  spacing  of  25  kHz.  Except 
for  a slight  degradation  due  to  Doppler  shifts  the  system  performance  for  data 
rates  of  2400  bps  and  frequency  spacings  of  3750  Hz  is  essentially  unchanged 
except  for  the  scaling  factor  in  time  and  frequency.  Since  the  effects  of 
Doppler  shifts  may  dominate  at  75  bps,  the  close  packing  of  user  signals  at 
the  teletype  rate  will  not  be  addressed  in  this  report. 

Typically  the  uplink  bandwidth  will  be  composed  of  a collection  of  sub- 
bands each  500  kHz  wide,  say,  containing  20  16  kbps  satellite  users  and  also 
some  interfering  LOS  users.  Supposing  that  the  LOS  users  produce  narrowband 
spurious  interference  they  are  modeled  as  narrowband-filtered  AWGN  at  25  kHz 
centers.  The  satellite  users  are  still  located  25  kHz  apart  but  in  between 
and  at  least  12.5  kHz  from  each  LOS  user. 

Because  of  the  natural  compactness  of  the  MSK  spectrum,  windows  were 
applied  at  the  transmitter  only  to  the  PSK  waveforms.  With  transmitter  win- 
dowing there  are  two  separate  cases  to  consider  depending  on  whether  the  trans- 
mitter is  or  is  not  peak-power -limited.  If  a non-rectangular  window  is  applied 
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and  the  transmitter  is  already  operating  at  full  power,  then  there  will  be  a 
signal  energy  loss  relative  to  the  channel  noise  at  the  receiver.  The  loss  is 
determined  by  the  ratio  of  the  energy  in  the  window  to  the  energy  in  a rectangu- 
lar window  of  the  same  length.  Simulated  results  for  QPSK  will  be  presented 
assuming  such  a peak-power  limitation,  and  just  for  one  particular  data  win- 
dow, for  which  this  loss  happens  to  be  3.07  dB.  That  is,  the  transmitter 
power  required  to  achieve  a certain  error  probability  and  E^/N^  with  QPSK  will 
be  up  to  about  3.1  dB  more  under  a peak- amplitude  restriction.  There  will  be  no 
such  loss  associated  with  the  constant  envelope  MSK  waveforms  since  a trans- 
mitter window  is  not  applied. 

Windows  were  applied  at  the  satellite  to  both  QPSK  and  MSK  waveforms.  In 
contrast  to  the  QPSK  case,  when  a satellite  window  was  applied  to  MSK,  it  was 

done  in  such  a way  that  no  intersymbol  interference  was  introduced. 

-3 

Symbol  error  probabilities  down  to  10  were  measured  by  computer  simula- 
tion for  uncoded  QPSK  at  16000  4-ary  symbols  per  sec  and  uncoded  MSK  at  both 
a 16  kbps  and  32  kbps  binary  symbol  rate.  The  rate  1/2  coded  performance  of 
QPSK  was  estimated  by  the  technique  explained  in  Section  4.5.  The  larger  MSK 
bit  rate,  corresponding  to  a rate  1/2  code,  proved  to  be  infeasible  for  the 
25  kHz  user  spacings,  so  coded  MSK  performance  was  not  estimated.  Uncoded 
MSK  at  the  16  kbps  rate  performed  quite  well,  however.  Thus,  the  performance 
of  uncoded  MSK  and  coded  QPSK  will  be  compared  on  the  common  basis  of  a 16  kbps 
information  rate  and  a 25  kHz  user  separation. 

The  required  E^/N^'s  for  a coherent  receiver  and  three  different  reliability 
levels  are  listed  in  Table  1.1.  The  parameter  E refers  to  the  received  power 
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TABLE  1.1 


SIMULATED  COHERENT  PERFORMANCE  WITH  WINDOWING  AND  ADJACENT  LOS  USERS  FOR  CODED  QPSK  AND  UNCODED  MSK 
a.  E = 10  dB 


Bit  Error  Probability 


E^/Nq  Required 


QPSK 


MSK 


Degradations  from  Theoretical 
for  AWGN 

QPSK 


MSK 


3-Bit  Decisions 

Hard  Decisions 

3-Bit  Decisions 

Hard  Decisions 

< 

10  2 

3.3  dB 

5.7  dB 

5.4 

dB 

0.8  dB 

1.6  dB 

0.2 

dB 

< 

O 

1 

U> 

4.5  dB 

7.2  dB 

9.2 

dB 

1.3  dB 

2.1  dB 

1.9 

dB 

< 

10-5 

6.8  dB 

9.0  dB 

_ 

2.0  dB 

2.4  dB 



b.  E = 15  dB 
Bit  Error  Probability 


E^/Nq  Required 


< 10 


-2 


< 10 


-3 


< 10 


-5 


QPSK 

3-bit  Decisions  Hard  Decisions 


4.2  dB 
6.1  dB 
11.1  dB 


8.6  dB 
12.3  dB 


MSK 

6.8  dB 
10.4  dB 


Degradations  from  Theoretical 
for  AWGN 

QPSK 

3-bit  Decisions  Hard  Decisions 


1.8  dB 

2.9  dB 
6.3  dB 


4.5  dB 
7.2  dB 


MSK 

1.6  dB 
3.1  dB 


Because  QPSK  requires  a transmitter  window,  up  to  3.1  dB  more  transmitter  power  is  needed  for  QPSK  than 
MSK  for  a given  E^/N^  if  the  transmitter  is  peak-power-limited. 


level  of  every  other  satellite  user  signal  relative  to  the  user  signal  being 
demodulated.  LOS  users  were  located  on  either  side  of  the  demodulated  satel- 
lite user  and  were  about  11  dB  above  a peak-power-limited  satellite  user  in 
pre-filtered  mean-square  signal  amplitude. 

Coded  QPSK  results  are  presented  for  both  3-bit  soft-decision  and  1-bit 
hard-decision  demodulation  in  the  satellite.  3-bit  decisions  have  an  advan- 
tage of  at  least  1/2  dB  and  2 dB  over  2-bit  soft-decision  and  hard-decision 
demodulation,  respectively,  for  the  convolutional  code  employed  [42],  With 
soft-decision  demodulation  in  the  satellite  and  end-to-end  coding  the  down- 
link bandwidth  necessary  will  generally  exceed  the  uplink  bandwidth  utilized. 
More  satellite  power  may  also  be  required.  The  additional  downlink  power  and 
bandwidth  required  would  depend  on  many  factors  including  the  modulation  for- 
mat, signaling  alphabet  size,  power  budget  margin,  etc.,  for  the  downlink  and 
therefore  will  not  be  treated  in  this  report.  If  it  turns  out  that  insufficient 
downlink  power  or  bandwidth  is  available  for  soft-decisions  and  end-to-end 
coding  then  either  hard-decision  demodulation  or  in-satellite  decoding  must 
be  used  with  QPSK. 

Observe  from  Table  1.1  that: 

CONSIDERING  THE  ADDED  3.1  dB  TRANSMITTER  POWER  REQUIRED  FOR 
WINDOWED  QPSK  UNDER  A STRICT  PEAK-POWER  LIMITATION  AT  THE 
TRANSMITTER,  MSK  IS  BETTER  THAN  CODED  QPSK  FOR  CVSD  MODULATED 
VOICE  (P,  < 10-2),  EVEN  FOR  SOFT-DECISION  DEMODULATION,  AND  FOR 

k _3 

VOCODED  SPEECH  (P,  < 10  ) FOR  HARD-DECISION  DEMODULATION, 

b 

MSK  has  the  additional  advantages  of  constant  envelope  signaling  and  not 
requiring  ISI  rejection  filters  and  coding  equipment.  Constant  envelope  sig- 
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naling  was  not  possible  for  QPSK  at  an  acceptable  performance  level  except 
with  near  perfect  power  control  (E  = 0 dB) . 

For  the  smaller  error  probabilities  of  Table  1.1: 

CODED  QPSK  WITH  SOFT  DECISIONS  IS  BETTER  THAN  MSK  FOR  VOCODED 

SPEECH  (P,  < 10"3)  OR  DATA  TRANSMISSIONS  (Pv  < 10_5)  . CODED 
k k _5 

QPSK  WITH  HARD  DECISIONS  SHOULD  OUTPERFORM  MSK  FOR  PL<  10 

b 

AND  E = 10  dB  INTERFERING  USERS,  EVEN  WITHOUT  A PEAK-POWER 
LIMITATION  AT  THE  TRANSMITTER. 

Although  accurate  performance  estimates  below  P^  = 10  are  lacking  for  MSK 
and  for  hard-decision  QPSK  at  E = 15  dB,  it  is  unlikely  that  MSK  could  out- 
perform QPSK  at  the  assumed  frequency  spacings  even  with  transmitter  window- 
ing and  higher  rate  codes.  Coded  QPSK  with  soft  decisions  would  be  substantially 
better  than  MSK  under  no  peak-power  limitation.  Unfortunately,  the  performance 
advantage  is  not  easily  estimated. 

Also,  from  the  performance  figures  in  Table  1.1: 

WITH  NO  PEAK-POWER  LIMITATION  AT  THE  TRANSMITTERS,  CODED  QPSK 

WITH  SOFT-DECISIONS  IS  SUPERIOR  TO  MSK  BY  AT  LEAST  2 dB  FOR 

P^  < 10-2,  AND  AT  LEAST  4 dB  FOR  PL  < 10_3,  WHILE  CODED  QPSK 
b b 2 

WITH  HARD  DECISIONS  IS  2 dB  BETTER  THAN  MSK  FOR  PL  < 10_J  AND 

b 

E = 10  dB  INTERFERING  USERS  BUT  WORSE  THAN  MSK  BY  ABOUT  2 dB 
FOR  E = 15  dB. 

In  summary,  then: 

WITH  THE  RELIABILITY  LEVELS  AND  PERFORMANCE  DEGRADATIONS  OF 
TABLE  1.1:  20  UNSYNCHRONIZED  16  kbps  SATELLITE  USERS  CAN  BE 

PACKED  INTO  A 500  kHz  UPLINK  BANDWIDTH  INCLUDING  NARROWBAND 
LOS  USERS  AND  COMMUNICATE  SIMULTANEOUSLY  IN  AN  FDMA  SATELLITE 
MODE  OF  OPERATION.  THE  RECEIVED  SIGNAL  STRENGTHS  MAY  VARY  BY 
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10  TO  15  dB  WITHOUT  THE  NEED  FOR  TRANSMITTER  POWER  CONTROL • 

BY  SCALING  THESE  RESULTS  UP  TO  7 UNSYNCHRONIZED  2400  bps  USERS 
SPACED  BY  3750  Hz  CAN  OCCUPY  A 25  kHz  BANDWIDTH  OR  AS  MANY  AS 
133  SUCH  USERS  CAN  COMMUNICATE  IN  A 500  kHz  BANDWIDTH. 

In  view  of  the  worst  case  crosstalk  analysis  for  the  close  packing  of 
MSK  waveforms  without  windowing  [8],  this  is  approximately  twice  the  number  of 
users  possible  in  these  bandwidths  without  the  special  processing  techniques 
developed  in  this  report. 

It  should  be  emphasized  that  computer  simulations  were  performed  with  a 
reasonable  but  quite  specific  model  for  LOS  interference  (see  Section  IV). 
However,  the  techniques  of  this  report  should  be  useful  in  dealing  with  other 
LOS  models.  For  example,  a longer  satellite  window  and/or  a sharper  cutoff 
frequency  could  be  used  against  stronger  and/or  wider  band  LOS  interference. 

The  net  increase  in  residual  intersymbol  interference  may  be  tolerable,  par- 
ticularly in  situations  where  transmitter  windowing  is  not  required.  This  may 
be  the  case  when  a single  satellite  user  is  allocated  a small  portion  of  a band- 
width occupied  predominantly  by  LOS  users. 
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II.  Analysis 

2.1  Mathematical  Model  for  an  MPSK  Digital  System 

In  this  section  a mathematical  model  for  the  signal  processing  parts  of  an 
M-ary  PSK  system  is  presented  and  some  terminology  is  established  which  will  be 
used  throughout  this  report.  Since  MSK  modulation  does  not  quite  fit  this 
model  and  is  not  considered  with  as  much  generality  as  MPSK  in  this  report,  the 
mathematical  model  for  MSK  will  be  treated  separately  in  Section  2.4. 

Figure  2.1  depicts  a digital  equivalent  model  for  an  MPSK  system  which 
omits  the  analog  and  hybrid  parts  of  a real  system  such  as  the  analog  to  digital 
converters,  analog  modulators,  etc.  Referring  to  that  figure  let  £ be  the  input 
symbol  index  for  a typical  user  u.  The  complex  sequence  of  input  symbols  to 
user  u,  XU£>  is  converted  into  a sequence  y^(n)  of  complex  numbers  by  using  a 
transmitter  data  window  t(n)  which  is  common  to  all  users.  As  an  example,  t(n) 
might  be  a sequence  of  constant  height  (referred  to  as  the  rectangular  window) 
which  spans  the  interval  between  the  arrival  of  x r and  x _ , - . This  particu- 
lar  window  has  the  distinct  advantage  of  permitting  constant  envelope  signaling 
for  a user  with  a peak-power-limited  transmitter,  but  on  the  other  hand  this  win- 
dow  does  little  spectral  shaping  and  thus  produces  potentially  serious  inter- 
user interference  or  crosstalk  (CT) . In  the  more  general  case,  y (n)  will  not 

u 

have  a constant  envelope  and  the  window  t(n)  will  span  more  than  one  symbol  in- 
terval in  order  to  reduce  the  spectral  splatter  of  user  u in  the  frequency 
bands  of  other  users.  However,  windows  that  span  more  than  one  symbol  interval 
introduce  intersymbol  interference  (ISI).  In  designing  a good  window,  the  min- 
imization of  CT,  the  minimization  of  ISI  and  producing  a signal  with  a nearly 
constant  envelope  are  conflicting  objectives.  This  will  become  clearer  in 
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Fig.  2.1. 


MPSK  digital  system  model  (only  one  user  path  is  shown) . 


what  follows. 


The  sequence  y (n)  is  modulated  by  the  sequence  , where  N is  the 

number  of  discrete  frequencies  in  the  digital  analysis  bandwidth  = exp(-27Tj/N) 
with  j = /-l,  and  is  an  integer  between  0 and  N — 1 corresponding  to  the 
carrier  frequency  for  user  u.  The  modulated  signal  is  then  converted  into  a 
continuous  signal  by  an  analog  to  digital  converter  (complex) , further  modulated 
by  an  uplink  carrier,  transmitted  over  a (presumed)  additive  white  Gaussian 
noise  (AWGN)  channel,  demodulated  to  an  intermediate  frequency  (IF)  and  re- 
sampled. This  whole  section  of  the  system  can  be  modeled  as  a digital  channel. 

In  this  channel  the  user's  signal  is  delayed  by  a delay  n^  with  respect 
to  the  other  users  (since  all  users  are  at  different  distances  to  the  satellite) 
and  attenuated  to  a real  amplitude  A^.  The  user's  frequency  is  also  perturbed 
by  an  amount  d^  due  to  Doppler  shifts  caused  by  relative  motion  between  the 
transmitter  and  satellite.  In  the  channel  the  perturbed  signal  is,  of  course, 
also  combined  with  similarly  modified  signals  of  other  satellite  users,  spuri- 
ous signals  from  line-of-sight  (LOS)  users  and  with  AWGN.  For  the  purpose  of 
this  model  and  the  simulation  of  Section  IV,  the  delay  and  Doppler  shift  of 
each  satellite  user  are  fixed  random  variables,  i.e. , n^  and  d^  are  randomly 
selected  at  the  start  of  an  experiment  but  held  constant  for  the  duration  of  a 
simulation  run. 

The  transmitted  signal  is  thus  perturbed  in  the  channel  to  form  the  modi- 
fied user  sequence 

(k  + d ) (n  — n) 

z (n)  = A y (n  — n )WM  u u u (2.1) 

u u u u N 

where  A^  is  a real  amplitude,  d^  is  a Doppler  shift  and  n^  is  a propagation 
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delay.  The  AWGN  sequence  will  be  denoted  by  w(n)  and  the  total  received  se- 
quence at  the  satellite  will  be  denoted  by  v(n) . 

At  the  satellite  the  signals  are  demodulated.  It  is  anticipated  that  this 

will  be  accomplished  in  U distinct  demodulation  channels  and  that  due  to  track- 
ing, the  center  frequency  k^,  Doppler  shift  d , and  delay  n^  for  each  user  sig- 
nal will  be  available  at  the  satellite.  In  the  special  case  that  the  user 

Doppler  shifts  and  delays  are  precorrected  so  that  the  effective  Doppler  shifts 
and  delays  are  zero,  then  the  demodulation  can  be  performed  on  all  of  the  chan- 
nels simultaneously . 

The  input  data  stream  at  the  satellite  is  segmented  into  subsequences  of 
length  by  means  of  a satellite  data  window,  s(n) . In  addition  the  satellite 
window  provides  further  spectral  shaping  to  help  reduce  crosstalk  and  reject 
LOS  interference.  This  window  is  applied  once  per  symbol  interval  on  each  de- 
modulation channel.  It  is  modified  to  compensate  for  the  Doppler  shift  d^  and 
it  is  delayed  so  that  it  is  applied  at  the  beginning  of  each  symbol  interval. 

If  f^(n)  denotes  the  output  sequence  from  the  windowing  operation  then 


d (n  — n ) 


fur,(n)  = s(n  “ nu>  v(n)WN 


, n < n < N + n -1 
U — — u 


As  a practical  matter  it  may  be  noted  that  the  multiplication  by  W 


. (2.2) 

d (n  — n ) 
u u 

N 

in  Eq.  (2.2)  need  be  performed  only  for  those  users  with  significant  Doppler 
shifts. 

* 

In  the  special  case  s(n)  = t (n)  the  satellite  and  transmitter  data  win- 
dows are  matched  and  the  satellite  window  then  functions  as  a matched  filter. 
Thus  in  this  case  the  signal-to-noise  ratio  at  the  output  of  the  DFT  is 
maximized. 
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A discrete  Fourier  transform  (DFT)  is  performed  once  each  symbol  interval 
for  each  user  to  separate  the  u-th  signal  from  the  others.  When  the  parameters 
of  the  DFT  are  properly  chosen,  one  obtains  as  the  k^-th  output  of  the  DFT  an 
estimate  of  the  input  symbol  x for  user  u and  input  symbol  interval 


N-l+n 


W = £ Wn)WN  (n"nu)ku  “]£  fur(n'  + nu}  Wn“  ^ ‘ (2‘3) 


N-l 


n'k 


n=n 


n =0 


Note  that  some  form  of  data  window  of  length  N is  implicitly  associated  with 

the  DFT  since  the  DFT  input  must  be  of  finite  length.  If  the  users  are  all  syn- 
chronized in  both  time  and  frequency  at  the  satellite,  s(n)  can  be  applied  just 

once  per  symbol  for  all  users  and  the  DFT's  needed  for  demodulation  can  all  be 
performed  simultaneously  using  a fast  Fourier  transform  (FFT)  algorithm  with 
some  savings  in  the  processing  hardware  (and  time)  required. 

The  output  of  the  DFT  may  or  may  not  represent  a good  approximation  to 
x depending  upon  the  windows  chosen  and  the  noise  level  in  the  channel.  Win- 
dows that  are  longer  than  the  intersymbol  distance  introduce  ISI  much  of  which 
can  be  removed  by  a simple  digital  filter  at  each  DFT  output.  Such  a filter. 


* 

as : 


Technically  this  is  non-standard  terminology.  The  DFT  is  normally  defined 


F(k) 


N-l 


f(n)wNnk 


, k = 0,  1, . . . , N-l 


In  general  it  converts  an  N-point  sequence  into  another  N-point  sequence.  In 
this  report  the  above  sum  is  often  computed  for  a single  value  of  k,  an  opera- 
tion which  will  be  called  a DFT  for  the  lack  of  a better  term.  The  fast 
Fourier  transform  (FFT)  algorithm  computes  this  sura  for  all  N values  of  k 
simultaneously. 
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with  impulse  response  h^,  produces  a complex  number  X^f  for  each  input  F^(k^) , 

A 

where  £ — £'  is  the  symbol  delay  of  the  filter.  Thus,  the  output  symbol  xu  ? 

is  an  estimate  of  the  earlier  transmitted  input  x . . 

u£ 

The  data  windows  t(n)  and  s(n)  are  restricted  to  be  finite  length  real 
discrete  time  functions  for  a less  complex  implementation  and  to  limit  the  ex- 
tent of  the  ISI.  It  is  assumed  that  the  transmitter  and  satellite  windows  ex- 
tend over  at  most  an  integral  number  r|  and  V of  symbol  intervals,  respectively, 
and  that  their  support  is  contiguous  and  centered  with  respect  to  those  inter- 
vals. These  assumptions  do  not  appear  to  limit  system  performance  significantly. 
Let  and  Ng  denote  the  length  of  t(n)  and  s(n),  respectively,  and  let  be 
the  intersymbol  distance.  ISI  will  be  introduced  whenever  either  window  is 

longer  than  the  intersymbol  distance  (N  > N or  N > N ) . 

t c s c 

One  way  to  conceive  of  the  operation  of  applying  the  transmitter  window 
is  to  first  generate  a sequence  x^(n)  from  the  x^  such  that 


x (n)  = 
u 


x _ , for  n = K N 
uC  c 

0,  otherwise 


(2.4) 


The  sequence  xu(n)  is  thus  a pulse  train  with  the  pulses  spaced  by  , the 
input  symbol  separation.  An  element  in  the  sequence  y^(n)  is  then  the  convolu- 
tion of  x^(n)  and  t(n).  This  procedure  is  illustrated  in  Fig.  2.2.  Mathemati- 
cally, y^(n)  can  be  written  as 

-E  <2-5> 

r 

All  but  n of  the  terms  in  the  sum  (2.5)  will  be  identically  zero  for  any  n due 
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xtl  (n) 


T<»-6-173i^ 


a.  A complex  signal  x^(n). 


tin) 


b.  A real  transmitter  window  t(n). 


Fig.  2.2.  Transmitter  windowing. 
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to  the  finite  length  of  t(n) . 


2.2  Effects  of  Transmitter  Windowing 

In  this  section  several  properties  of  the  transmitter  window  are  discussed. 

First,  it  is  noted  that  for  equimagnitude  symbols,  |xu^|  = say,  con- 
stant envelope  signaling  is  possible  only  for  a rectangular  transmitter  window 
that  spans  exactly  one  symbol  interval.  A rectangular  window  of  greater  sup- 
port results  in  non-constant  envelope  transmission,  as  can  be  deduced  from 
Eq.  (2.5). 

The  nN^-point  DFT  of  an  N^-point  rectangular  transmitter  data  window  of 
height  h is  given  by 


T(k)  = h 


exp 


/ tt(N  -1)  \ 

(-J  k) 


, k = O',  1, . . . , r)N  — 1 
c 

(2.6) 


A graph  of  | T ( k) | reveals  that  by  modifying  the  value  of  Nt  the  magnitude  and 
width  of  the  main  lobe  of  the  function  can  be  altered  but  little  control  over 
the  sidelobes  is  available.  Hence,  only  limited  control  over  the  interference 
among  users  is  possible  with  this  window.  Furthermore,  making  the  window 
longer  by  increasing  beyond  increases  ISI.  Using  a more  general  shaped 
window  one  can  improve  both  the  spectral  shaping  due  to  the  window  for  reducing 
CT,  and  the  ISI,  but  this  must  be  done  at  the  expense  of  the  signal  energy  for 
a peak-power-limited  transmitter.  Since  the  signal-to-noise  ratio  is  reduced 
with  a shaped  window,  the  probability  of  error  due  to  the  AWGN  in  the  channel 
increases.  There  are  thus  four  sources  of  error:  line-of-sight  (LOS)  inter- 

ference, Gaussian  channel  noise  (GN) , crosstalk  (CT),  and  intersymbol 
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interference  ( IS I) . By  varying  window  shapes  and  lengths,  one  can  find  various 
operating  points  which  trade  off  the  relative  mix  of  these  sorts  of  errors. 

When  choosing  windows,  however,  all  four  sources  of  error  must  be  kept  in  mind. 

Consider  next  the  average  signal  energy  penalty  for  a peak  power  limited 
transmitter  in  AWGN  with  a non-rectangular  window.  Assuming  statistically 
independent  zero-mean  input  symbols  the  total  mean  received  signal  energy  from 
user  u during  one  symbol  interval  is 


E 

uc 


where 


(2.7) 


N -1 
s 

<}>  (£)  s(m)  t(m  - AN  ) l = 0,  ± 1,  ± 2, . . . (2.8) 

m=0 

is  the  crosscorrelation  function  between  the  two  windows,  and  the  shift  is 
measured  in  symbols. 

Clearly,  for  fixed  A , N , N and  for  a transmitter  window  which  is  peak- 

u t c 

amplitude -limited  with  amplitude  h,  the  window  yielding  the  largest  mean  energy 
is  rectangular  with  constant  amplitude  h.  The  mean  received  energy  from  user 
u for  such  a window  is 

E = N (A  h)2  (2.9) 

u t u 

2 

by  Eqs.  (2.7)  and  (2.8).  It  is  convenient  to  choose  h such  that  E = A , i.e., 

u u 

so  that  for  the  rectangular  window  of  support  N^,  $^(0)  = 1*  With  any  other 
window,  the  energy  received  from  user  u will  be  less.  In  fact  with  this  nor- 
malization the  average  fraction  of  energy  received  from  user  u with  respect  to 
E^  is  simply  Euc/Eu  = 4>tt:(0)  £ 1 from  Eqs.  (2.8)  and  (2.9).  Thought  of  another 
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way,  there  is  at  most  an 


Lt  - - 10  log10<J>tt(0)  dB  (2.10) 

loss  in  mean  received  signal  energy  due  to  normalization  at  a peak-power-limited 
transmitter.  For  most  of  the  non-rectangular  windows  tried  in  the  simulation 
L was  in  the  range  1.89  to  3.07  dB . 

2.3  Combatting  Intersymbol  Interference 

There  are  a number  of  ways  to  handle  ISI,  and  there  is  a vast  literature 
dealing  with  time-varying  dispersive  channels.  However,  because  of  the  rather 
special  nature  of  the  ISI  which  appears  here,  a simple  linear  recursive  filter 
of  low  order  is  used  to  reduce  ISI  [23,43,44].  The  approach  is  suboptimal,  but 
its  performance  in  this  context  compares  favorably  with  the  optimal  approach 
using  the  Viterbi  algorithm,  and  it  is  much  simpler  to  implement.  Because  this 
scheme  was  so  simple  and  yet  performed  nearly  optimally,  other  schemes  were 

not  tried  since  there  seemed  to  be  little  to  gain  from  alternate  approaches 
[23,45,46]. 


A reasonable  and  tractable  criterion  for  designing  the  digital  filter  of 
Fig.  2.1  is  that  of  minimizing  the  mean- square  error  (MSE) 


e = X — x 
u£  u£  u£' 


(2.11) 


between  the  complex  filter  output  symbol  X for  user  u and  the  corresponding 
system  input  symbol  x^.  I n fact,  one  is  really  interested  in  minimizing  the 
decision  error  probability  for  a given  signal-to-noise  ratio.  But  this  is  a 
much  more  difficult  if  not  intractable  analytical  task.  Another  approach  might 
be  to  design  the  digital  filter  to  maximize  the  signal-to-noise  ratio  (SNR)  at 
the  filter  output.  This  latter  approach  has  been  tried  but  for  typical 
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interference  models  (delays,  Doppler  shifts  and  amplitudes)  the  optimum  SNR  and 
minimum  MSE  (MMSE)  filters  are  not  significantly  different. 


Although  the  digital  filter  design  will  be  based  on  the  MMSE  criterion, 
system  performance  will  be  measured  by  computer  simulation  to  estimate  deci- 
sion error  probabilities.  A comparison  will  be  made  between  the  suboptimal 
MMSE  filter  and  the  optimal  Viterbi  algorithm  approaches  for  reducing  ISI.  Al- 
though theoretical  receiver  operating  characteristics  are  known  for  the  Viterbi 
algorithm  in  AWGN,  the  formulas  are  not  easily  applied  so  performance  will  also 
be  measured  by  simulation. 

2.3.1  MMSE  Recursive  Filter  Design 
The  MMSE  digital  filter  that  operates  once  per  channel  symbol  to  reduce 
ISI  is  designed  as  follows.  For  simplicity  the  subscript  ufs  are  temporarily 
deleted.  if  F denotes  the  input  to  the  filter  and  X denotes  its  output,  the 
filter  is  defined  by  the  difference  equation 


N 


N 


X 


e ai  fn  , + C - i ^ 

i=0  d 


i=l 


i XC  - i 


(2.12) 


F^  comes  from  Eq.  (2.3)  and  includes  the  effects  of  crosstalk  (CT) . The  filter 
design  problem  is  the  problem  of  finding  the  {a^,  B_^}.  In  general  these 
parameters  can  be  complex  but  if  the  data  windows  are  real,  the  {a^,  B^}  must 
be  real.  > 0 is  the  filter  symbol  delay.  and  are  the  number  of  zeros 

and  poles,  respectively,  of  the  digital  filter. 

Differentiation  of  Eq.  (2.11)  with  respect  to  each  a and  B^  using  Eq.  (2.12) 
and  setting  each  partial  derivative  to  zero  yields  the  nonlinear  system  of  equa- 
tions 
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^ai(RFF(i“£)  + RFF(£_i))  +?ei(RXF(£"Nd_i)  + RFX(Nd+i‘£)) 

= (Rxp(£-Nd)  + RFx(Nd-£)),  A-0,l,---,Nz  (2.13a) 

Z^d-m)  + Rxx(m_i))  +?ai(RXF(i-Nd_m)  + RFX(Nd+m-i)) 


(RxX^m^  + ’ m = 1 >2,  * * * ,N 


(2.13b) 


where  R.  (j-k)  = p.a  [47]  for  the  complex  random  variables  p and  a.  These 
po  j k 

equations  are  nonlinear  because  the  correlation  functions  R^(£)  > R^p(£)  9 
RpX(£)  , ^xX(^)  and  Rxx(£)  are  t^iemselves  aH  functions  of  the  {o^,  8^}.  The 
fact  that  the  windows  and  the  (a^,  8^}  are  real  guarantees  that  all  of  the  cor- 
relation functions  in  Eq.  (2.13)  are  real  functions  of  £.  Since  for  any  cor- 
relation  function  R^(£)  = R (-£)  the  equations  of  (2.13)  can  be  simplified 
to  the  following: 


^aiRFF(Jl-i)  + S6.RFX(Nd+  i-£)  = RFx(Nd-A),  A-0,1,---, 
i i 

( 

?aiRFX^Nd+  m"i^  +X^BiRxx(m-i)  =RKx(“m)»  m- 1,2,  • • • ,N  . 


N 

z 

(2.14) 


It  should  be  noted  that  the  matrix  of  coefficients  in  Eq.  (2.14)  is  symmetric, 
considering  the  {a^,  8^}  as  unknowns. 

The  nonlinear  equations  of  (2.14)  were  solved  iteratively.  First  the 
system  of  linear  equations 
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(2.15) 


a.RFFa-i)  = Rpx(Nd  - £),  l = 0,  1,...,  Nz 
i 

was  solved  to  give  an  initial  vector  } . Using  the 

digital  filter  impulse  response  (a  , = ())  new  correlation  functions 

were  then  computed  and  the  full  set  of  Eqs.  (2.14)  was  solved  for  and 

The  filter  hj^(ot^\  J3^)  was  then  used  to  compute  a new  set  of  coefficients 
(2)  (2) 

or  ' and  J3V  } , and  so  on.  Excellent  convergence  of  the  {a^,  3^}  has  been  ob- 
served within  just  a few  iterations  (typically  3 to  10)  . 

* 

Filters  were  computed  for  several  values  of  N , and  and 

and  13  dB , where  E,  is  the  mean  received  energy  per  bit  from  user  u and  N is 
b o 

the  power  spectral  density  of  the  complex  zero-mean  AWGN  process.  Good  results 
were  usually  obtained  for  = 2 to  4 and  = 0 or  1.  Only  very  slight  im- 
provement was  observed  for  larger  values  of  or  (consistent  with  [23]). 

Thus,  very  simple  linear  filters  are  adequate  for  mitigating  most  of  the 
ISI,  crosstalk  (CT) , and  GN  for  the  windows  tested.  A canonical  filter  config- 
uration with  two  zeros  and  one  pole  is  shown  in  Fig.  2.3,  where  F(z)  and  X(z) 
are  taken  as  the  z-transforms  of  F and  X , respectively.  Recall  that  although 
the  input,  output  and  contents  of  the  delay  elements  are  complex,  the  multi- 
pliers are  real.  It  should  probably  be  emphasized  at  this  point  that  if  the 
length  of  each  window  is  less  than  the  intersymbol  distance,  then  the  ISI  is 
zero  and  no  digital  filters  are  necessary,  since  their  primary  function  is  to 
remove  ISI.  There  will  be  ISI  if  either  window  spans  more  than  one  symbol. 

At  this  point  it  is  helpful  to  write  a more  explicit  expression  for  Fr 
using  Eqs.  (2.1)  - (2.5)  and  Fig.  2.1: 


It  is  assumed  that  E^/N^  and  L^_  are  both  expressed  in  dB. 
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F(z) 


H(z) 


aQ  + a(/z  + a2/z 
1 -^/z 


Fig.  2.3.  Canonical  2-zero,  1-pole  linear  digital  filter. 
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Fr(ku>  = ^WN^ku + du^ns(n)  (w(n  + iO 
^ n 


(2.16) 


+ SAiWN(ki  + di>  (Vn-nU)Exi,c+  r t(n+nu-n.-rNc))  . 


The  various  components  of  Eq.  (2.16)  will  now  be  isolated.  As  mentioned  before, 
in  addition  to  the  signal  there  are  three  types  of  noise  excluding  EOS  inter- 
ference — interuser  noise,  intersymbol  interference  and  channel  noise.  The 
signal  is  that  part  of  Eq.  (2.16)  which  is  due  solely  to  the  C-th  input  symbol 
from  user  u: 


fSIG  * VSt<0)  V 


(2.17) 


The  channel  interference  (GN)  results  from  the  processed  channel  noise  w(n) : 


GN  =2w(n+nu)s(n)WN(ku+du)n> 


n 


(2.18) 


The  intersymbol  interference  is  due  to  the  other  x 


u?* 


frcr  = A 2^1  X rJ_  <i>  _(r). 
^SI  U r?0  u>£+r  st 


(2.19) 


The  interuser  noise  is  due  to  everything  else: 


fCT  = f Ai  ^Xi^+r  ^S(n)  t(n+nu-ni"rNc) 
lfu  r n 


(2.20) 


(ku+du)n  + (k±+d . ) (n  -n-n  ) . 
N 
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Thus,  Fr(k  ) = f + f + f + fnrr  for  the  £-th  symbol  interval  of  user  u. 
G U bib  GN  ibl  Cl 

From  Eqs.  (2.17)  - (2.20)  one  may  calculate  the  following  autocorrelation 
functions  for  the  signal  and  noise  terms: 


rsig(^  ^ <Vst(0»  RSIG^0^ 


(2.21) 


under  the  assumption  that  R (^  — CT)  = 

XX  £ G ( U , 


e = V 

otherwise’ 


RGN(£)  = Vss^  = Vss(£)  ’ 


(2.22) 


risi(£) 


Au  4r*st(r)*st(r  + A) 
r?0 

r^-i 


(2.23) 


and 


rct(*>  = ^2  Ai  ^2  ^2  t^n  + nu  ~ ni  - rNc^ 

i^u  r n,m 

* t(m  + n - n - (r  + £)N  ) WM  (2.24) 

u i c N 


assuming  statistically  identical  but  independent  users. 

In  contrast  to  the  results  Eqs.  (2.21)  - (2.23),  R (£)  is  complex  in 

c i 

general.  Sufficient  conditions  for  a real  R (£)  are  equiamplitude  users  sym- 

c 1 

metrically  spaced  in  frequency  with  respect  to  user  u,  i.e.,  the  k^'s  are  sym- 
metric about  k^  modulo  N,  with  zero  delays  and  Doppler  shifts.  These  condi- 
tions are  not  unreasonable  at  least  on  the  average  and  therefore  the  filter  is 

designed  to  these  conditions.  A real  R ( £)  in  the  design  simplifies  both  the 

c x 

solution  of  the  nonlinear  Eq.  (2.14).  and  the  resulting  filter  by  permitting 


half  the  number  of  real  coefficients  (ol,  8^1-  This  follows  from  Eqs.  (2.12) 
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and  (2.14)  and  the  fact  that 


rffU)  = rsigw  + rgn(£)  + Risia)  + rct(£)  * (2.25) 

In  this  report  the  digital  filter  is  non-adaptive  because  the  {a^,  8^} 
coefficients  do  not  change  with  time. 

A couple  of  issues  involving  the  recursive  filter  that  have  been  con- 
sidered in  some  detail  will  be  mentioned  only  briefly  in  the  remainder  of  this 
section. 

The  first  issue  has  to  do  with  prefiltering  at  the  transmitters  for  re- 
ducing ISI.  So  far  the  MMSE  digital  filters  have  been  located  in  the  satellite, 
i.e.,  at  the  output  of  each  user  demodulation  channel.  Conceivably,  these 
compensating  filters  could  instead  be  placed  at  the  input  of  each  user  terminal. 
The  design  problem  for  this  input  filter  option  is  very  similar  to  that  just 
presented  for  the  output  filters,  although  the  equations  are  slightly  differ- 
ent from  Eq.  (2.13). 

After  working  out  the  necessary  details  for  an  input  filter  design,  one 
is  convinced  of  the  feasibility  of  the  approach  but  reluctant  to  program  the 
scheme  for  computer  simulation.  The  potential  gains  from  an  input  filter  sim- 
ulation do  not  seem  to  justify  such  an  investigation  at  present.  Furthermore, 
constant  envelope  signaling  at  the  transmitter  is  impossible  with  the  input 
filter. 

The  second  issue  has  to  do  with  estimating  the  irreducible  error  probabil- 
ity due  to  crosstalk  and  ISI.  The  signal-to-noise  ratio  (SNR)  for  the  overall 
system  can  be  expressed  in  terms  of  the  autocorrelation  functions  of  the  var- 
ious signal  and  noise  components  and  the  transfer  function  of  the  MMSE  filter. 
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Since  the  symbol  error  probability,  P^,  in  the  presence  of  AWGN  can  be  calcu- 
lated from  the  SNR,  the  theoretical  P^  for  a typical  simulation  was  estimated 
using  just  the  SNR.  Unfortunately,  this  estimate  was  quite  pessimistic  be- 
cause the  noise  is  neither  white  nor  Gaussian  and  because  delays  and  Doppler 
shifts  can  corrupt  the  overall  system  performance.  As  a consequence  extensive 
simulation  appears  to  be  necessary  to  establish  an  accurate  performance  floor 
for  large  SNRfs. 

2.3.2  The  Viterbi  Algorithm  Approach 

The  optimal  but  more  complex  maximum  likelihood  decoding  approach  to  the 
ISI  problem  is  now  considered. 

When  the  transmitter  and  satellite  windows  are  identical,  the  basic  sys- 
tem configuration  of  Fig.  2.1  conveniently  fits  the  Forney  model  [48]  for  the 
mitigation  of  ISI  using  the  Viterbi  algorithm  (VA) . In  order  to  greatly  sim- 
plify the  model  it  is  assumed  that  the  interuser  interference  can  be  treated 
as  AWGN.  Simulation  results  indicate  that  this  assumption  is  not  a bad  one 
particularly  since  windows  can  be  selected  to  make  crosstalk  (CT)  negligible. 

From  Eqs.  (2.17)  - (2.19)  the  received  symbol  sequence  at  the  output  of 
the  DFT  of  Fig.  2.1  can  be  written  as 

f(r)  = fSIG(r)  + fisi(r)  + fGN(r)  = 4>ts(r)  + fGN(r)  (2.26) 

where  © represents  discrete  convolution.  In  the  z-transform  domain  Eq.  (2.26) 
becomes 

F(z)  = fsig(z)  + FISI(z)  + FGN(Z)  = M(z)  $ts(z)  + FGN(z)  • (2'27) 

The  Gaussian  noise  f_.,(r)  has  the  autocorrelation  function  N cj)  (r)  from 

GN  o ss 
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Eq . (2.22),  and  the  power  spectrum  N $ (z).  Since  4>  is  an  autocorrelation 

OSS  s s 

function  its  z-transform  can  be  factored  as 

<S>  (z)  = P ( z)  P(l/z)  (2.28) 

ss 

where  all  the  zeros  of  P(z)  are  inside  the  unit  disc  of  the  z-plane.  Since 

(J)  is  a sequence  of  length  2v  — 1,  P(z)  has  V — 1 zeros  and  is  the  z-transform 
s s 

of  a sequence  of  length  V,  p(r) . In  the  case  that  the  transmitter  and  satellite 
windows  are  identical,  t(n)  = s(n)  and  Eq . (2.27)  can  be  rewritten  as 

F(z)  = AX(z)  P(z)  P(l/z)  + N(z)  P ( 1 / z ) (2.29a) 

W(z)  = AX(z)  P(z)  + N(z)  (2.29b) 

where  n(r)  is  an  AWGN  sequence.  Equation  (2.29b)  was  obtained  from  Eq.  (2.29a) 
by  division  by  P(l/z).  In  terms  of  an  implementation  w(r)  is  the  output  of  a 
digital  filter  with  impulse  response  p ^(r)  when  f(r)  is  the  input.  Since  the 
noise  component  of  Eq.  (2.29b)  is  white  this  filter  will  be  referred  to  as  a 
whitening  filter. 

The  next  step  is  to  attempt  the  removal  of  the  ISI  factor  P(z)  in  Eq. 
(2.29b)  and  thereby  obtain  the  input  symbol  sequence  with  the  transform  X(z). 

The  VA  performs  a maximal  likelihood  estimation  of  the  entire  input  sequence 
in  AWGN  and  is  quite  well  suited  to  combatting  ISI.  With  this  scheme  the  out- 
put digital  filter  and  decision  rule  of  Fig.  2.1  is  replaced  by  the  whitening 
filter  and  the  VA. 

The  question  of  the  necessary  degree  v — 1 of  the  P polynomials  is  con- 
sidered next.  Simulations  suggest  that  matched  transmitter  and  satellite  win- 
dows that  extend  over  more  than  two  symbols  may  not  be  required  except  when  the 
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interfering  users  are  very  strong  and/or  very  closely  packed  in  frequency 
relative  to  the  satellite  user  being  demodulated.  For  the  more  reasonable 
parameter  values  of  the  simulation,  the  additional  CT  protection  was  not  deemed 
worth  the  increase  in  ISI  and,  in  the  case  of  the  peak-power  limitation,  the 
signal  energy  penalty  with  respect  to  AWGN.  Because  of  this  and  the  fact  that 
there  is  no  ISI  unless  the  windows  span  more  than  one  symbol  interval,  a speci- 
fic instance  of  the  VA  is  described  and  simulated  only  for  the  case  V = 2. 
Although  the  VA  implementation  would  become  more  complex,  there  are  no  problems 
with  extending  the  method  to  larger  values  of  V. 

For  V = 2,  Eq.  (2.28)  has  the  form 

$ (z)  = a/z  + b + az  = (c/z  + d)  (cz  + d)  (2.30) 

ss 

where  a = (p  (1)  = <J>  (-1),  b = <J)  (0)  and  P(z)  is  the  monomial  in  1/z.  The 

s s s s s s 

whitening  filter  has  the  form  shown  in  Fig.  2.4a.  In  the  absence  of  noise  n(r) 
the  output  W(z)  of  the  whitening  filter  of  Fig.  2.4a  equals  the  output  of  the 
feedforward  filter  of  Fig.  2.4b  with  input  AX(z)  and  response  P(z). 

In  the  VA  all  possible  outputs  (y(r)} 

y (r)  = d • x(r)  + cx(r  — 1)  (2.31) 

of  the  noiseless  feedforward  filter  of  Fig.  2.4b  at  a given  symbol  time  r are 
computed  for  all  possible  current  input  symbols  (x(r) } and  all  possible  im- 
mediately preceding  input  symbols  {x(r  — 1)}.  These  computations  need  only  be 
done  once,  so  the  answers  could  be  precomputed  and  stored  in  a read-only 
memory  (ROM).  A subsequence  of  these  (y(r)}  is  then  compared  to  the  actual 
received  sequence  z(r)  out  of  the  whitening  filter  in  such  a way  as  to  decide. 
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W(z) 


F (z) 
P(1/z) 


a. 


+ ► Y(z)  =AX(z)P(z) 


H(z)  = P(z)  = c/z  + d 


b. 

Fig.  2.4.  Whitening  filter  (a)  and  Viterbi  algorithm  feedforward  filter  (b). 
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with  some  fixed  symbol  delay,  on  the  most  likely  input  sequence  x(r).  Naturally, 
the  most  likely  input  sequence  would  produce  the  y(r)  subsequence  that  is 
closest  to  the  received  sequence  z(r)  in  an  Euclidean  distance  sense. 

The  VA  is  usually  described  [49]  in  terms  of  a trellis  diagram  and  a set 
of  state  sequence  paths  of  finite  length  L as  depicted  in  Fig.  2.5  for  M = 4. 

For  each  symbol 'time  r there  are  MV  ^ possible  states  or  nodes  in  the  r-th 
column  of  the  trellis,  one  node  for  each  of  the  possible  input  symbols  (x(r)}. 

Let  G^(r)  be  a real  number  associated  with  the  node  in  column  r corresponding 
to  the  input  symbol  x^(r).  For  the  first  value  of  r,  G^(r)  = 0 for  each  value 
of  i.  Let  i1  be  the  row  index  for  the  nodes  in  column  r — 1.  Assume  that  there 
are  M surviving  paths  of  length  L,  where  the  head  of  each  path  is  a distinct 
node  in  column  r — 1 and  where  the  tail  node  of  each  path  is  a (not  necessarily 
distinct)node  in  column  r — L.  The  likelihood  of  the  path  headed  by  the  sym- 
bol x^,(r  — 1)  is  measured  by  the  number  G^,(r  — 1).  The  first  step  of  the 
VA  is  to  determine  the  latest  transition  to  column  r for  each  path  in  a new 
set  of  M surviving  paths  as  follows.  The  tail  nodes  of  the  current  set  of 
surviving  paths  are  deleted.  Then  for  each  i the  M transition  numbers  {g^,^) 


is  to  the  received  z(r),  the  smaller  the  g^,^.  The  latest  transition  for  the 


8i'i  = * xi^r^  + cxi<(r  ~ D)|2 


(2.32) 


are  computed  using  Eq.  (2.31),  one  g^,^  for  each  possible  transition  from  a 
node  in  column  r — 1.  The  closer  the  computed  y(r)  = d • x^(r)  + cx^,(r  — 1) 


new  surviving  path  headed  by  node  ir  is  determined  by  an  index  i'  satisfying 


Git (r)  = Gi(r) , where 


(2.33) 
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The  other  M-l  transitions  to  node  ir  are  discarded  as  indicated  by  the  dashed 
lines  to  node  2r  in  Fig.  2.5.  The  next  step  of  the  VA  is  to  output  the  input 
symbol  corresponding  to  the  tail  node  of  the  new  most  likely  surviving  path. 
The  head  of  this  path  is  specified  by  the  index  i satisfying  G(r)  = G^(r), 
where 


_ , N min  ^ , v 
G(r)  = . G.(r) 


(2.34) 


The  decoded  symbol  is  associated  with  the  node  in  column  r-L-1  belonging  to 
the  most  likely  path.  Thus,  there  is  a symbol  delay  of  L-l  due  to  decoding. 

In  general,  the  complexity  of  the  VA  varies  as  MV  where  M is  the  number 
of  possible  input  symbol  values  and  V is  the  length  of  the  impulse  response 
of  the  feedforward  filter  with  transform  P(z) . No  multiplications  are  required 
during  the  operation  of  the  VA  to  obtain  the  y(r)'s,  for  a generalization  of 
Eq.  (2.31)  can  be  precomputed  and  stored  in  a memory  of  size  MV  2q-bit  words, 
where  2q  is  the  number  of  bits  required  to  store  a complex  number.  MV  complex 
subtractions,  2MV  real  multiplications  and  MV  real  additions  are  required  per 
symbol  interval  to  compute  the  {g^,^},  however,  cf.  Eq.  (2.32).  Another 
M^(V  ^ r^al  additions  and  MV  ^ (MV  ^)-way  comparisons  are  required  to  compute 
the  real  numbers  (G^(r)}  assigned  to  the  heads  of  the  new  paths,  cf . Eq.  (2.33), 
and  one  (M^^-way  comparison  is  needed  for  Eq.  (2.34).  Furthermore,  assuming 
that  M is  a power  of  two,  LMV  ^ (v-1)  log^  M-bit  words  of  memory  are  required 
to  store  the  state  sequence  paths  of  length  L. 

The  total  costs  of  the  VA  with  complex  arithmetic  are  indicated  in  Table 
2.1.  As  already  mentioned  V=2  is  the  primary  value  of  interest.  For  quadri- 
phase  shift  keyed  (QPSK)  modulation,  a principal  candidate  for  the  communica- 


tions system  at  hand,  the  x(r)  symbol  alphabet  is 


± 1 * 3 


/2 


>,  so  M=4.  A path 
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TABLE  2.1 


COMPLEXITY  OF  THE  VITERBI  ALGORITHM 


Operation 


Number  of  Operations 
per  Symbol 


real  q-bit  multiplications 
real  q-bit  additions 
real  (M^  ^)-way  comparisons 


2MV 

MV_1  + 1 


+ 3M 


.V 


Memory 

V 

M 2q-bit  words 

LMV  ^ ( V— 1)  log^  M-bit  words 


length  L of  10  to  20  is  typical  for  good  performance.  A q of  5 is  usually  ade- 
quate, though  a 4 or  8 bit  quantization  is  more  convenient;  a word  length  that 
is  a multiple  of  four  bits  is  preferable  for  ease  of  implemention.  For  these 
parameter  values,  then,  the  VA  requires  32  real  multiplications,  64  real  addi- 
tions and  5 4-ary  comparisons  per  symbol.  A ROM  of  16  8-or  16-bit  words  and 
a random-access  memory  (RAM)  of  4 20- to  40-bit  words  is  required. 

In  comparing  the  complexity  of  the  VA  and  MMSE  approaches  to  the  ISI 
problem,  the  whitening  filter  of  v-1  delay  elements,  cf.  Fig.  2.4a,  must  be 
included  in  the  cost  of  the  VA  approach,  of  course.  Considering  the  fact  that 
a typical  MMSE  linear  recursive  filter,  cf.  Fig.  2.3,  requires  8 real  multi- 
plications and  6 real  additions  per  symbol,  and  a RAM  of  only  3 8-or  16-bit 
words,  it  is  obvious  that  the  VA  approach  is  much  more  complex  than  the  MMSE 
approach.  In  the  computer  simulations  the  performance  of  these  two  schemes 
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was  quite  comparable  (see  Section  IV). 

2.4  Mathematical  Model  for  an  MSK  Digital  System 

As  stated  at  the  beginning  of  Section  II,  minimum  shift  keying  (MSK)  is 
not  treated  with  as  much  generality  as  MPSK  in  this  report.  In  particular,  no 
special  windowing  is  applied  at  the  transmitters  with  MSK  and  the  satellite 
windows  are  restricted  to  those  which  introduce  no  ISI.  Under  these  conditions 
there  is  no  need  for  a compensating  digital  filter  or  a Viterbi  algorithm  ap- 
proach for  reducing  ISI.  These  restrictions  on  the  windows  have  yielded  pro- 
mising simulated  results,  nevertheless,  because  of  the  relatively  small  cross- 
talk (CT)  between  unsynchronized  and  unwindowed  MSK  waveforms  [8]. 

Figure  2.6  depicts  a digital  equivalent  model  for  an  MSK  system  with  much 
the  same  notation  as  used  in  Fig.  2.1  for  an  MPSK  system.  In  the  MSK  case  two 
sequences  of  + 1 input  symbols  to  user  u,  x^  anc*  xus^>  are  convert:ec*  into  a 
sequence  y^(n)  °f  complex  numbers 


/ \ TTn  . TTn 

y (n)  = x cos  -r—  + j x sin 
u ucC  2N  j us4  2N 

c c 


(2.35) 


where  2N^  is  the  intersymbol  distance  in  each  input  sequence.  The  values  of 
xuc^  are  permitted  to  change  only  when  n = N (modulo  2N^) , while  the  values 
of  x may  change  only  when  n = 0 (modulo  2N^) . 

Although  the  expressions  for  z^(n)  and  f ^ (n)  given  by  Eqs.  (2.1)  and 
(2.2)  still  apply,  MSK  is  demodulated  differently  than  MPSK.  Referring  to 


Fig.  2.6,  the  outputs  F ^(k  ) and  F (k  ) are  computed  from 
° c^  u sC  u 
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.EN 


Fig.  2.6. 


MSK  digital  system  model 


(only  one  user  path  is 


shown) . 


(2.36a) 


N -1 
c 


^ nk 

F (k  ) = > f r(n  + n ) Wxt  U 

cC  u u£  u N 


COS 


TTn 

2N 


n=-N 


2N  -1 
c 


nk 


‘sC 


(k  ) = f (n  + n ) W U sin 

u UC  u N 


n=0 


u • 7111 

2N 


(2.36  b) 


where  the  range  of  the  summation  index  for  F (k  ) and  F (k  ) corresponds  to 

c£  u sc;  u 

the  interval  when  x „ and  x ^ must  be  constant,  respectively.  Symbol  estimates 

ucC  usC 

are  formed  as 


ucC 


+1,  if  Re{F  (k  )}  > 0 
c£  u 

-1,  otherwise 


(2.37a) 


us£ 


+1,  if  Im{F  (k  ) } > 0 

SC,  u 


(2.37b) 


-1,  otherwise 

The  various  components  of  Eq.  (2.36b)  for  the  sine  demodulation  channel 
will  now  be  isolated.  A similar  exercise  could  be  performed  for  the  cosine 
demodulation  channel.  Using  Eqs.  (2.1)  and  (2.2)  and  excluding  LOS  users  in 
Fig.  2.6,  the  sine  channel  signal  due  solely  to  the  £-th  sine  input  symbol 


x _ from  user  u is 
usC 


CsSIG  Xus£  Au 


2N  -1  0 

c 2 

r \ r • 7111  \ 

s (n)  (sm  “2jj~) 

n=0 


(2.38) 


The  satellite  window  s(n)  is  limited  implicitly  to  the  intersymbol  distance 

2N  for  MSK,  so  there  is  no  ISI.  The  sine  channel  interference  due  to  the 
c 


noise  w(n)  is 


2N  -1 

( c (k  +d  )n j 

fsGN  ’ Im  j E ”<n  + nu>  S(n)  (Sin  lb)  UN  “ j 


(2.39) 
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and  the  sine  channel  crosstalk  is  given  by 


2N  -1 


cf.  Eqs . (2.18)  and  (2.20),  respectively. 


W, 


(k  +d  ) n + (k . + d ) (n . -n-n 
u u l i i u 


(2.40) 
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III.  WINDOWS 


The  use  of  data  windows  for  bandwidth  conservation  in  an  FDMA  mode  is  the 
central  theme  of  the  approach  taken  in  this  report.  This  section  is  devoted  to 
the  choice  and  design  of  the  data  windows  chosen  for  the  system  environment  out- 
lined in  Section  I. 

3.1  System  Aspects 

For  purposes  of  exposition  and  for  the  bulk  of  the  simulations,  an  infor- 
mation rate  of  16  kbps  was  assumed.  It  was  further  assumed  that  PSK  modulation 
would  be  used,  that  the  center  frequencies  of  adjacent  channels  would  be  25  kHz 
apart  and  that  a maximum  contiguous  uplink  bandwidth  of  approximately  500  kHz 
is  available  for  satellite  users.  Thus  up  to  20  users  can  be  accommodated. 

A rate  1/2  convolutional  code  for  greater  reliability  and  QPSK  modulation  are 
presumed,  although  the  following  discussion  also  applies  to  an  uncoded  BPSK 
system.  In  either  case  the  channel  symbol  rate  is  16000  symbols/sec  which 
corresponds  to  an  intersymbol  spacing  T^  of  62.5  \\sec.  Recognizing  that  nar- 
rowband LOS  users  may  also  have  frequency  allocations  in  the  uplink  bandwidth, 
the  worst  case  situation  of  20  LOS  users  spaced  by  25  kHz  and  interspersed  be- 
tween 20  satellite  users  is  presumed.  The  problem  is  thus  to  design  data  win- 
dows to  control  both  the  interuser  and  LOS  interference  with  respect  to  any 
satellite  user. 

A discrete  window  sequence  is  used  to  achieve  spectral  shaping,  both  to 
minimize  interference  between  adjacent  satellite  users  and  to  reduce  inter- 
ference from  LOS  users  operating  at  nearby  frequencies.  In  the  worst  case  the 
LOS  users  will  be  operating  on  25  kHz  centers  and  the  satellite  users  will  also 
be  forced  to  operate  on  25  kHz  centers  in  the  cracks  between  LOS  channels  as 
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shown  in  Fig.  3.1.  Clearly,  the  available  center  frequency  to  cutoff  bandwidth 
for  a given  satellite  user  must  not  exceed  12.5  kHz.  In  order  to  see  what 
effect  the  choice  of  a transmitter  window  t(n)  has  upon  the  shape  of  the  signal 
spectrum,  the  signal  at  the  output  of  one  particular  user  transmitter  is 
analyzed. 

Let  £^(n)  represent  the  discrete  sequence  formed  by  the  application  of  the 
transmitter  window.  £^(n)  can  then  be  written,  cf.  Fig.  2.2  and  Eq.  (2.5),  as 

Ri 

?u(n)  = E xur  exp[j  ¥ kun]  t(n  “ rV  * (3-1) 

r"Ro 

Here  and  are  used  to  denote  indices  of  the  first  and  last  symbols  trans- 
mitted, respectively,  and  N is  a parameter  which  controls  (along  with  k^)  the 
center  frequency  of  the  satellite  user.  If  the  window  t(n)  is  of  finite  length 
(the  only  case  considered),  then  only  a few  of  the  terms  in  Eq.  (3.1)  contri- 
bute to  the  sum  for  a given  value  of  n.  The  Fourier  transform  of  this  signal 
can  be  written  as 

00 

- / j(i)x  ~ / \ -i^n 

- (e  ) = > £ (n)e  J 

u ' JLj  ^u 

n=-°° 


"2  E xur  exp|j  ¥ kunJ  t(n_  rNc)  exP[-j^n] 

n=-oo  r=RQ  L 


2iTk 

1 

/ 27Tk  v-i 

|exp  |j 

V 

1 

3 

)Ev exp 

l'r-R0 

[-JrNc  (“  — )] 

(3.2) 

where  TCe^)  is  the  Fourier  transform  of  t(n).  The  power  density  spectrum  of 
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Fig.  3.1.  Worst  case  frequency  allocations. 
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this  signal  can  be  expressed  as 


|5(ejW)  |2  = |T(exp  [ j (w 


(3.3) 


assuming  that  the  information  symbols  are  uncorrelated. 

The  D/A  converter  of  Fig.  1.1  is  clocked  with  an  intersample  time  of 
Tg  = T^/N^.  Thus,  if  H(jf2)  denotes  the  transform  of  the  continuous  (analog) 
output  of  the  D/A  converter,  then  from  the  sampling  theorem 

jftT 

z (jfl)  = Ts  5 (e  s)  (3.4a) 

and 

|H  (jfi)|2  = T2  |T[exp  j (OTs  . (3.4b) 


This  is  an  important  result  for  it  identifies  the  shape  of  the  signal 

spectrum  with  the  Fourier  transform  of  the  transmitter  window.  In  particular 

the  power  spectral  density  of  the  u-th  user  is  the  magnitude  of  the  Fourier 

transform  of  the  window  squared  and  displaced  in  frequency  by  k /NT  . Thus, 

to  locate  satellite  users  at  odd  multiples  of  12.5  kHz  we  require 

1/N  = 12.5  x 10^  T . Since  T = T /N  with  T =62.5  ysec,  this  implies  that 
s see  c 


N = 


25N 


c 32 

For  the  number  of  samples  per  symbol  to  be  an  integer,  N should  thus  be  a 
multiple  of  32.  N represents  the  total  number  of  satellite  channels  available, 
but  only  half  of  these  channels  can  have  center  frequencies  which  are  located 
at  odd  multiples  of  12.5  kHz.  In  order  to  have  at  least  20  usable  16  kbps 
channels  and  to  permit  a simple  FFT  realization  in  the  satellite,  N was  chosen 
to  be  the  least  power  of  two  larger  than  32,  i.e.,  64.  This  fixes  the  number 
of  samples  per  symbol  at  50  and  yields  a total  analysis  bandwidth  of  800  kHz 
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which  represents  the  basic  clock  rate  for  the  A/D  and  D/A  converters.  A simi- 
lar analysis  could  be  carried  out  for  different  data  rate  channels. 

Two  different  classes  of  windows  were  selected  for  study  — the  Dolph- 
Chebyshev  windows  and  the  prolate  spheroidal  windows.  For  a given  length  and 
cutoff  frequency  each  of  these  windows  is  optimum  in  some  sense  — a Dolph- 
Chebyshev  window  has  a minimum  equiripple  response  outside  of  the  main  lobe  of 
its  frequency  response  and  a prolate  spheroidal  window  has  minimum  energy  out 
of  band.  Each  of  these  windows  has  an  adjustable  parameter  so  that  its  band- 
width can  be  controlled,  and  each  has  a real,  even  frequency  response  which 
aids  in  the  implementation  of  the  system.  Since  the  prolate  spheroidal  window 
was  both  easy  to  work  with  and  also  optimal  in  a particularly  useful  way,  other 
windows  were  not  used  extensively.  It  is  extremely  doubtful  that  another  win- 
dow can  be  found  that  will  improve  the  performance  of  the  system  measurably. 

These  two  window  families  will  be  discussed  in  the  following  two  sections. 

Since  a given  window  could  be  applied  at  either  the  transmitter  or  the  satellite, 
the  generic  notation  w(n)  for  the  data  window  and  WCe^00)  for  its  Fourier  trans- 
form will  be  used  just  in  these  sections. 

3.2  The  Dolph-Chebyshev  (DC)  Window 

The  DC  window  represents  the  optimal  equiripple  window.  A DC  window  of 

length  N with  cutoff  frequency  0)  has  a digital  frequency  response 
w c 

cos(N  cos  1 (z  cos((jo/2))) 

M/  10)N  W o r x 

W(eJ  ) = — (3.5a) 

cosh  (N  cosh  z ) 
w o' 

z = 7 j-r-  (3.5b) 

o cos  (a)  /2) 

c 

as  shown  in  Fig.  3.2.  For  a given  sidelobe  height  6,  the  DC  window  has  the 
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Fig.  3.2. 


Dolph-Chebyshev  spectral  window. 


smallest  raainlobe  width  co^  of  any  window  of  the  same  order.  Alternatively, 
for  a given  mainlobe  width  the  DC  window  has  the  smallest  maximum  sidelobe  height. 
All  of  the  sidelobes  are  of  the  same  height,  so  the  window  is  called  equiripple. 
This  implies  that  to  a given  user  all  other  equal  power  users  contribute  the 
same  amount  of  interference  whether  they  are  located  on  an  adjacent  channel  or 
separated  by  several  channels.  The  total  crosstalk  energy  received  from  other 
users  is  thus  directly  proportional  to  the  number  of  such  users. 

For  a communication  channel  with  a center  frequency  to  cutoff  bandwidth 
of  f^  = 00^/271  Hz  and  an  intersaraple  time  Tg,  according  to  Helms  [50]  the  side- 
lobe  ripple  6 is  bounded  by 


log 


2 

10  6 


f T N 
< C s w 

0.73 


(3.6) 


Using  T =1.25  ysec  for  a given  f , one  can  compute  the  tradeoff  between  the 
s c 

window  length  and  the  peak  sidelobe  ripple.  This  information  is  displayed 
graphically  in  Fig.  3.3. 

In  order  to  adequately  reduce  interference  among  different  users  for  the 
signal  strengths  expected,  and  system  parameters  described  above,  more  than  20  dB 
of  rejection  will  be  required.  Consequently,  DC  windows  will  have  to  be  more 
than  50  samples  (one  symbol  length)  long  as  seen  from  Fig.  3.3.  This  will  re- 
sult in  intersymbol  interference  which  will  have  to  be  mitigated  by  post  fil- 
tering in  the  satellite. 

Several  Dolph-Chebyshev  data  windows  and  their  frequency  responses  were 
computed  and  utilized  in  the  simulations.  It  was  noted  that  short  windows 
possess  a large  spike  in  their  time-domain.  This  is  a possibly  severe  draw- 
back if  short  windows  were  to  be  used  with  peak-power-limited  transmitters  for 
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Fig.  3.3.  Sidelobe  level  vs  window  length  for  Dolph-Chebyshev  window  and 
analysis  bandwidth  of  800  kHz. 
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it  limits  the  average  power  which  can  be  transmitted. 

Each  data  window  was  computed  by  periodically  sampling  the  frequency  re- 
sponse of  Eq.  (3.5a)  and  then  computing  the  inverse  discrete  Fourier  trans- 
form (DFT) 


w(n) 


(3.7) 


The  data  windows  were  normalized  for  a peak-power-limited  transmitter.  In 

particular  they  were  normalized  so  that  the  energy  in  a rectangular  window  of 

length  N is  unity,  cf,  Eqs.  (2.8)  -(2.10).  A non-rectangular  window  represents 
w 

an  energy  loss  due  to  the  peak-power-limitation  at  the  transmitter.  In 
those  cases  where  there  is  no  peak-power-limitation,  this  energy  loss  does  not 
appear. 

The  sampled  crosscorrelation  functions  (p  (£) 

st 

VN  -1 
c 

<p  *.(£)  = s(m)  t (m  — £N  ) 

st  / ^ c 

for  some  of  the  computed  DC  data  windows  are  displayed  in  Fig.  3.4.  The  center 
pulse  (£  = 0)  from  each  set  represents  the  fraction  of  the  signal  energy  which 
remains  after  the  two  windowing  operations.  The  height  of  the  pulses  on  the 
sides  represents  the  amount  of  intersymbol  interference  (ISI)  which  is  present. 
Only  when  both  windows  span  only  one  symbol  is  no  ISI  introduced.  Figure  3.4a 
depicts  the  case  with  windowing  at  both  the  transmitter  and  the  satellite. 

Here  s(n)  = t(n)  = w(n)  using  Eqs.  (3.5)  and  (3.7).  Figure  3.4b  and  c de- 
pict DC  windowing  only  at  the  satellite  with  a rectangular  shaped  data  window 
at  the  transmitter.  In  Fig.  3.4a  and  b the  transmitter  and  satellite  windows 
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Fig.  3.4.  Crosscorrelation  functions  for  Dolph-Chebyshev  windows  (transmitter 
window  t(n)  of  length  and  satellite  window  of  length  Ng) . 
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are  aligned  when  £ = 0 , i.e.,  they  have  a common  vertical  centerline.  In 
Fig.  3.4c  the  non-zero  portion  of  the  transmitter  window  is  shifted  one  symbol 
interval  to  the  right.  The  cases  represented  in  Fig.  3.4b  and  Fig.  3.4c  permit 
constant  envelope  signaling  at  the  transmitter,  because  t(n)  is  rectangular  and 
one  symbol  in  length,  while  allowing  the  option  of  centering  the  non-zero  part 
of  a satellite  window  of  length  c Ng  2N^  over  a 3 or  2 symbol  interval, 
respectively. 

3.3  The  Prolate  Spheroidal  (PS)  Window 

The  prolate  spheroidal  windows  represent  optimal  windows  in  the  sense  that 
of  all  windows  with  length  they  have  the  smallest  amount  of  energy  outside 
the  band  (-0)  , ^ ) • For  this  reason  they  are  superior  to  the  Dolph- 
Chebyshev  (DC)  family  in  many  cases.  They  are  a little  more  difficult  to  gen- 
erate, however.  Their  generation  will  be  more  completely  understood  if  the 
form  of  the  window  is  derived  from  basic  considerations. 

Consider  any  window  w(n)  of  length  with  Fourier  transform 

N -1 

W(ejw)  = ^ w(n)  e^Wn  . (3.8) 

n=0 

One  wishes  to  find  that  w(n)  for  which  the  energy  contained  in  the  frequency 

band  (—03  , u)  ) is  a maximum, 
c c 

Restated,  the  problem  is  to  find  the  w(n)  which  maximizes 
0) 

rc  . 2 

J |w  (eJW)  | do) 

6 = — . (3.9) 

" _ 2 
J |W  (ejw)|  du>- 
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This  problem  is  similar  to  the  problem  for  continuous  signals  which  was  solved 
by  Landau  and  Pollack  using  prolate  spheroidal  wave  functions  [51,52].  The 
following  approach  for  discrete  signals  obviously  borrows  from  their  earlier 
work. 

From  Eq.  (3.8)  and  Parseval's  relation 
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(3.10) 


Thus,  Eq.  (3.9)  becomes 
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where  w (m)  is  the  complex  conjugate  of  w(m).  From  the  Cauchy-Schwartz 
inequality, 
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with  equality  if  and  only  if  b(n)  = c * a(n),  where  c is  any  constant.  If 
Eq.  (3.12)  is  applied  to  Eq . (3.11)  with 


a(n)  = w(n) 


N -1 
w 

1 T — ' * 

b(n)  = — 23  w 
m=0 


sin  a)  (m-n) 
c 

(m-n) 


then  the  maximum  value  of 


will  occur  when 


N -1 
w 


V w(n)  = \ ' w 
m=U 


(m) 


sin  0)  (n-m) 

c 


TT(n-m) 


(3.13a) 


(3.13b) 


(3.14) 


Here  A^  is  a constant,  and  w(m)  is  real.  According  to  the  form  of  Eq.  (3.14) 

2 

Aq  must  thus  be  an  eigenvalue  of  that  equation.  The  value  of  3 that  results 

is  straightforwardly  computed  as  A^  by  substituting  Eq.  (3.14)  into  Eq.  (3.11). 

2 

Thus  to  maximize  3 one  must  find  the  largest  eigenvalue  of  Eq.  (3.14)  and  choose 
w(n)  to  be  the  corresponding  eigenfunction. 

Some  examples  of  these  PS  windows  obtained  by  a standard  numerical  method 
[53]  are  shown  in  Fig.  3.5.  The  data  windows  are  normalized  in  the  same  way 
as  the  DC  data  windows.  The  energy  loss  L^_  for  peak-power-limited  transmitters 
is  also  shown  in  Fig.  3.6. 

In  Fig.  3.7  the  fraction  of  out-of-band  energy  1-A  is  plotted  as  a func- 
tion of  N for  different  values  of  user  bandwidth.  In  contrast  to  the  DC  win- 
w 

dows  the  co-channel  interference  decreases  as  the  channel  separation  is  increased 
due  to  the  spectral  roll-off  of  the  PS  windows.  Thus,  with  PS  transmitter  win- 
dows system  performance  will  be  less  ' sensitive  to  interference  from  equal 
power  users  at  larger  frequency  deviations. 
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a.  50-point  data  window. 


Fig.  3.5.  Prolate  spheroidal  windows  for  f =12.5  kHz  and  an  analysis 
bandwidth  of  800  kHz. 


59 


c.  76-point  data  window. 


Fig.  3.5  Continued. 
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Fig.  3.5  Continued. 
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h.  128-point  spectral  window. 
Fig.  3.5  Continued. 
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Fig.  3.6.  Maximum  energy  loss  for  prolate  spheroidal  windows  with  peak- 
power-limited  transmitter . 


Fig.  3.7.  Fraction  of  out-of-band  energy  for  prolate  spheroidal  windows. 
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The  crosscorrelations  <f>  (£,)  for  these  PS  data  windows  are  displayed  in 

st 

Fig.  3.8.  This  figure  is  the  prolate  spheroidal  counterpart  to  Fig.  3.4. 

Detailed  descriptions  of  how  these  windows  performed  will  be  deferred  to 
Section  IV,  although,  as  might  be  expected,  the  PS  windows  are  to  be  preferred 
over  the  DC  windows.  The  PS  windows  are  more  natural  for  use  in  conjunction 
with  the  MMSE  filters  of  Section  2.3.  In  addition,  the  decreasing  sidelobe 
level  of  the  PS  family  is  probably  of  greater  importance  than  the  slightly 
narrower  mainlobe  of  the  DC  family.  Finally,  the  spikes  in  the  DC  data  window 
can  be  a disadvantage  for  peak-power-limited  transmitters. 
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IV.  PERFORMANCE 


In  this  section  selected  results  from  an  extensive  computer  simulation 
using  data  windows  for  bandwidth  conservation  are  presented. 

4.1  Rationale  of  Simulation 

Simulation  was  used  to  estimate  symbol  error  probabilities  since  analytic 
expressions  for  the  error  are  difficult  to  obtain  in  the  presence  of  interuser 
interference.  Even  when  the  windows  shape  the  user  spectra  sufficiently  so 
that  the  crosstalk  (CT)  between  users  is  negligible  and  computable  bounds  for 
the  intersymbol  interference  (ISI)  error  can  be  applied,  simulation  has  still 
proven  useful  to  determine  how  much  ISI  can  be  removed  with  a simple  MMSE  fil- 
ter. Direct  simulation  has  the  additional  advantage  that  all  signaling  schemes 
can  be  compared  on  a common  basis. 

A flexible  set  of  computer  programs  has  been  developed  for  generating 
Dolph-Chebyshev  (DC)  and  prolate  spheroidal  (PS)  windows,  designing  MMSE  filters 
for  given  transmitter  and  satellite  windows,  and  performing  the  indicated  simu- 
lations for  satellite  users  employing  BPSK,  QPSK,  and  CPFSK  modulation.  The 
scheme  of  Forney,  cf.  Section  2.3.2,  which  uses  the  Viterbi  algorithm  to  obtain 
a maximum  likelihood  estimate  of  the  transmitted  data  sequence  in  the  presence 
of  AWGN  and  ISI  has  also  been  programmed.  LOS  users,  when  simulated,  are 
modeled  as  colored  Gaussian  noise  generated  by  passing  AWGN  through  a narrow- 
band  digital  filter. 

Although  many  options  can  be  simulated,  computer  time  is  expensive.  A 
typical  experimental  run  over  1000  symbols  consumes  on  the  order  of  1 minute 
of  IBM  370/168  CPU  time.  Such  run  lengths  are  sufficient,  however,  to  permit 
reasonably  tight  confidence  intervals  for  the  measured  error  probabilities  in 
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the  range  of  2 to  12  dB.  The  authors  have  chosen  to  spend  computer  time 

investigating  a variety  of  system  options  with  moderate  statistical  variances 
in  the  resulting  data  instead  of  just  a few  options  with  longer  runs  and  smal- 
ler variances. 

4.2  Qualitative  Results 

Initially  LOS  users  were  ignored  and  data  were  gather  for  QPSK  satellite 

users  under  several  different  windowing  modes.  Allowing  for  the  case  that  some 

user  signals  could  be  as  much  as  10  dB  weaker  than  others  at  the  receiver,  it 

was  quickly  established  that  non-rectangular  windows  spanning  more  than  one 

symbol  were  necessary  for  a viable  level  of  performance  (e.g.,  a P symbol  error 

-2  -3 

probability  between  10  and  10  at  an  E_/N  of  10  dB)  . Hence,  with  QPSK  users 

b o 

it  was  not  possible  to  achieve  the  system  objectives  without  introducing  ISI. 
Next,  it  was  determined  that  constant  envelope  (CE)  QPSK  signaling  (rectangular 
transmitter  window  spanning  one  symbol  interval)  and  a long  satellite  window 
covering  more  than  one  symbol  (which  introduces  ISI)  performed  rather  poorly. 

The  superior  CT  properties  of  MSK  suggest  that  CE  MSK  signaling  with  either  a 
rectangular  or  non-rectangular  satellite  window  s(n)  may  be  feasible.  In  fact 
simulations  indicate  that  MSK  modulation  is  a very  attractive  keying  scheme  for 
uncoded  data.  Finally  matched  PS  windows  (s(n)  = t(n))  yielded  promising  re- 
sults for  QPSK  with  window  lengths  between  one  and  two  symbol  intervals.  In 
this  latter  case  a simple  MMSE  digital  filter  was  as  effective  at  removing  in- 
tersymbol interference  as  a maximum  likelihood  estimator  based  on  the  Viterbi 
algorithm.  The  DC  windows  proved  to  be  uniformly  inferior  to  the  corresponding 
PS  windows. 

Although  BPSK  (binary  phase  shift  keying)  users  were  not  simulated 
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extensively,  the  existing  data  and  theory  suggests  that  the  qualitative  obser- 
vations from  the  study  of  QPSK  would  also  apply  to  BPSK.  The  major  difference 
is  that  the  BPSK  simulation  reflects  the  uncoded  situation  while  the  16000 
4-ary  symbols  per  sec  rate  QPSK  allows  for  a performance  improvement  via  a rate 
1/2  convolutional  code.  (See  Section  4.5.) 

In  the  absence  of  interfering  LOS  users  the  windows  for  the  satellite  users 
were  somewhat  arbitrarily  designed  to  have  a cutoff  frequency  of  12.5  kHz.  It 
was  natural  to  consider  a smaller  cutoff  frequency  in  dealing  with  LOS  inter- 
ference. This  provides  an  extra  degree  of  freedom  that  hopefully  would  help 
reduce  LOS  interference  by  narrowing  the  main  lobe  of  the  spectral  window. 

For  a fixed  window  length,  however,  a narrower  main  lobe  means  larger  sidelobes 
hence  increased  CT  among  users.  The  simulations  showed  that  the  original  cut- 
off frequency  of  12.5  kHz  was  superior  even  in  the  presence  of  the  LOS  inter- 
ference modeled,  evidently  because  the  sidelobe  effects  dominated. 

Before  presenting  some  quantitative  results  to  illustrate  the  foregoing 
conclusions,  the  format  of  the  simulations  will  be  described. 

4.3  Scope  of  the  Simulation 

Twenty-one  prolate  spheroidal  and  Dolph-Chebyshev  windows  have  been  gen- 
erated for  simulation.  The  window  characteristics  are  summarized  in  Appendix  A. 
Window  lengths  between  50  and  200  samples  were  used  spanning  between  one  and 
four  input  symbols. 

Windows  whose  lengths  were  between  64  and  100  samples  were  simulated  in 
two  ways  — once  with  the  total  span  extending  over  only  two  symbol  intervals, 
and  once  with  the  span  extending  over  parts  of  three  symbol  intervals.  The 
centering  of  a window  is  not  important  in  the  design  of  a MMSE  compensating 
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filter  but  the  simulations  for  the  two  cases  were  different.  Although  the 
simulation  evidence  is  not  conclusive,  the  authors  tend  to  favor  using  three 
symbol  intervals  for  these  windows. 

Simulations  were  also  performed  using  50-point  windows  which  due  to  their 
one  symbol  length  produce  no  intersymbol  interference.  For  these  windows  no 
correction  filters  are  necessary  and  the  output  of  the  DFT  is  itself  an  esti- 
mate of  the  transmitted  symbol. 

The  cutoff  frequency  for  the  spectrum  of  the  overall  signal  was  varied 
between  8 and  12.5  kHz.  The  windows  with  lower  cutoff  frequencies  were  used 
in  experiments  with  simulated  LOS  users.  In  the  bulk  of  the  simulations, 
however,  LOS  users  were  not  simulated. 

It  was  always  assumed  that  several  users  would  be  transmitting  over  sep- 
arate channels  simultaneously.  Due  to  a variety  of  factors  (atmospheric  con- 
ditions, unequal  distances  from  the  transmitters  to  the  satellite,  etc.)  these 
users  will  not  be  received  at  the  satellite  with  equal  energies.  It  is  impor- 
tant nonetheless  that  the  weaker  channels  still  be  reliable.  Thus,  it  was 
generally  assumed  that  all  users  except  the  one  being  demodulated  would  be  re- 
ceived with  energy  E times  that  of  the  user  of  interest.  E varied  from  0 to 
20  dB.  On  some  runs  the  E for  a given  user  was  taken  to  be  a random  variable 
uniformly  distributed  between  0.25  and  2.25.  Interfering  users  were  either 
chosen  to  be  synchronized  or  unsynchronized.  In  the  latter  case  each  user  was 
given  a random  delay.  These  delays  were  uniformly  distributed  between  0 and  1 
symbol  interval.  They  were  also  given  Doppler  shifts  relative  to  one  another 
of  up  to  + 500  Hz.  In  all  cases  the  user  signal  being  demodulated  is  assumed 
to  be  synchronized  with  respect  to  its  own  receiver. 
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The  characteristics  of  the  MMSE  filters  for  reducing  ISI  are  listed  in 
Appendix  B.  The  coefficients  of  these  filters  depend  upon  the  specific  trans- 
mitter and  satellite  windows  chosen.  Each  filter  was  designed  assuming  10  equal 
energy  satellite  users,  zero  delays  and  Doppler  shifts,  and  a certain  level  of 
channel  noise. 

The  satellite  user  frequencies  employed  in  the  filter  design  and  simula- 
tion were  either  closely  packed  (which  represents  a worst  case  crosstalk  situa- 
tion) or  loosely  spaced.  In  the  packed  case  adjacent  users  were  placed  at 
frequencies  + 25  kHz,  + 50  kHz,...  from  the  user  being  demodulated.  In  the 
spaced  case,  the  frequency  spacing  was  three  times  as  great. 

Although  all  of  the  digital  filters  were  designed  to  minimize  the  mean- 
square  error  (MSE)  between  the  actual  and  desired  outputs,  the  filters  could 
contain  any  number  of  poles  and  zeros.  The  MSE  figure  was  used  to  select  pro- 
mising filter  configurations.  The  smaller  the  MSE,  the  greater  the  degree  of 
ISI  removed.  This  figure  generally  increases  with  the  channel  noise  and  with 
the  window  length  as  shown  in  Fig.  4.1  for  some  2-zero,  1-pole  filters  de- 
signed for  matched  PS  windows.  Adding  more  than  one  pole  (feedback  path)  to 
the  filter  made  little  significant  improvement  in  either  the  MSE  or  in  the 
performance  of  the  system.  Adding  more  than  two  zeros  (three  feedforward 
terms)  to  the  filter  reduced  the  MSE  significantly  only  for  windows  longer 
than  100  samples. 

4.4  Quantitative  Data 

The  most  significant  simulation  data  justifying  the  qualitative  results 
of  Section  4.2  will  be  presented  in  this  section.  In  almost  all  cases  each 
simulated  datum  was  based  on  1000  observations. 
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Fig.  4.1.  MSE  f igure-of-merit  for  filter  designs  with  peak-power-limited 
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In  Fig.  4.2  the  performance  of  a 50-point  prolate  spheroidal  window  at 
the  satellite  is  shown  both  with  a constant  envelope  (CE)  50-point  rectangular 
transmitter  window  and  with  a 50-point  matched  PS  transmitter  window.  When 
the  users  are  all  synchronized,  the  matched  window  is  clearly  superior  and  in 
the  unsynchronized  case  neither  window  performs  adequately.  These  data  were 
gathered  for  10  dB  interfering  users  (E  = 10)  and  were  plotted  assuming  a 
strictly  peak-power-limited  transmitter  (loss  = 1.96  dB  due  to  transmitter 
window).  The  effects  of  Doppler  shifts  were  unsignificant  in  this  case. 

There  is  an  advantage  to  using  a CE  transmitter  window  against  channel 
noise,  when  the  transmitter  is  peak-power-limited.  In  Fig.  4.3,  therefore,  a 
50-point  CE  transmitter  window  is  paired  with  a 128-point  PS  satellite  window 
and  a 150-point  DC  satellite  window.  The  other  users  in  this  case  had  random 
timing  asynchronies,  but  no  significant  Doppler  shifts.  The  cases  of  0 dB, 

6 dB  and  10  dB  interfering  users  are  shown.  It  will  be  observed  that  the  PS 
window  is  slightly  better  than  the  DC  in  this  case  and  that  neither  window  is 
adequate  in  the  presence  of  strong  interfering  users.  These  results  are  sup- 
ported by  other  simulations.  From  these  data  it  can  be  safely  concluded  that 
nearly  perfect  power  control  would  be  required  for  CE  signaling  with  QPSK. 

On  the  other  hand  there  was  a clear  improvement  for  unsynchronized  users  with 
the  longer  windows  over  the  case  presented  in  Fig.  4.2. 

Neither  the  one  symbol  matched  window  configuration  nor  the  CE  transmitter 
window  just  discussed  will  perform  adequately  if  the  synchronization  and  trans- 
mitter power  control  requirements  are  relaxed.  By  contrast,  for  QPSK,  matched 
satellite  and  transmitter  windows  which  span  more  than  one  symbol  interval  will 
perform  well  in  this  situation  as  evidenced  by  the  data  presented  in  Fig.  4.4. 
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Fig.  4.2.  Uncoded  16k  symbols/sec  QPSK  performance  for  one-symbol  prolate 
spheroidal  satellite  window. 
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Fig.  4.3.  Uncoded  16k  symbols/sec  QPSK  performance  for  constant  envelope 
transmitter  window  and  long  satellite  window. 
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Uncoded  16k  symbols/sec  QPSK  performance  for  matched  prolate 
windows  with  peak-power-limited  transmitters. 
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For  equal  amplitude  satellite  users  the  data  agree  very  closely  with  the 
theoretical  QPSK  curve  reflecting  the  maximum  3.07  dB  energy  penalty  in  the 
presence  of  AWGN  due  to  the  non-rectangular  signal  envelope.  The  compensating 
MMSE  digital  filter  performs  well  in  reducing  the  ISI  introduced  by  these 
relatively  short  76-point  (about  3/2  symbol  intervals)  windows.  These  short 
windows  do  a good  job  in  controlling  crosstalk  since  there  is  only  a 0.5  to 
2.5  dB  degradation  from  the  theoretical  curve  for  interfering  users  as  much  as 
10  dB  stronger  than  the  user  being  demodulated. 

The  use  of  a whitening  filter  and  Viterbi  algorithm  to  remove  intersymbol 
interference  yielded  eventually  the  same  performance  as  the  MMSE  filter  for  an 
E^/Nq  up  to  about  12  dB . The  former  maximum  likelihood  approach  is  to  be  pre- 
ferred for  higher  signal-to-noise  ratios. 

With  longer  windows  there  is  increased  ISI  and  an  increased  penalty  to  be 
paid  for  the  shape  of  the  signal  envelope.  With  200-point  matched  PS  windows 
the  maximum  energy  penalty  is  3.76  dB  and  there  is  a 1 dB  increase  in  the 
residual  ISI  in  the  E^/No  range  of  6 to  12  dB.  Thus,  in  general,  as  the  data 
windows  are  made  longer  the  CT  is  reduced  at  the  price  of  increased  ISI  and 
vulnerability  to  AGWN  under  a peak-power  constraint  at  the  transmitters. 

Figure  4.5  illustrates  the  consistent  observation  that  matched  DC  windows 
are  inferior  to  matched  PS  windows,  especially  when  the  interfering  satellite 
users  are  strong  relative  to  the  user  being  demodulated.  There  is  also  more 
residual  ISI  with  a DC  window  than  with  the  corresponding  PS  window.  This 
discrepancy  is  particularly  severe  for  long  windows.  Consequently,  Dolph- 
Chebyshev  windows  will  be  eliminated  from  further  consideration. 

Up  to  now  only  satellite  users  have  been  simulated.  Next  some  results 
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Fig.  4.5.  Uncoded  16k  symbols/sec  QPSK  performance  for  matched  Dolph- 
Chebyshev  windows  with  peak-power-limited  transmitters. 
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obtained  with  LOS  as  well  as  satellite  users  will  be  presented. 

The  designers  of  the  satellite  system  presumably  have  no  real  control  over 

the  LOS  users  but  only  limited  knowledge  of  the  LOS  frequency  allocations,  nomi- 

* 

nal  bandwidth  and  modulation  formats,  etc.  A typical  LOS  user  would  be  of 
relatively  narrow  bandwidth  and  located  at  a multiple  of  25  kHz.  Aside  from 
this  a given  LOS  modulation  characteristics  could  be  just  about  anything. 
Accordingly,  the  LOS  users  are  modeled  simply  as  narrowband-filtered  (colored) 
Gaussian  noise  at  the  appropriate  center  frequencies.  The  nominal  bandwidth 
of  these  users  was  set  in  the  range  3 to  5 kHz  and  their  pre-filtered  mean- 
square  signal  amplitude  was  set  about  8 dB  and  11  dB  above  the  user  of  interest 
for  QPSK  and  MSK,  respectively , assuming  a rectangular  transmitter  window.  In 
these  simulations  two  LOS  users  were  assumed  located  on  each  side  of  the  user 
being  demodulated  as  shown  in  Fig.  4.6.  Although  more  LOS  users  could  be  in- 
cluded in  this  model,  the  filtering  action  of  the  satellite  window  guarantees 
that  any  additional  LOS  users  would  have  only  a secondary  effect. 

Figure  4.7  presents  data  for  three  different  pairs  of  matched  windows. 
Window  2 has  a frequency  cutoff  of  12.5  kHz,  window  15  cuts  off  at  8 kHz  and 
window  18  cuts  off  at  10  kHz.  All  three  of  the  PS  windows  are  76  points  long. 

It  should  be  noted  that  the  12.5  kHz  window  performs  better  despite  the  fact 
that  it  has  the  largest  maximum  energy  penalty  L^_ . In  comparing  Figs.  4.4  and 
4.7  it  is  seen  that  the  presence  of  these  LOS  users  does  not  degrade  perfor- 
mance significantly  for  in  the  range  4 to  12  dB  for  this  window. 


There  is  a predominance  of  amplitude  modulated  voice  communication  with 
a 6 kHz  bandwidth  at  UHF. 
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Fig.  4.6.  Line-of-sight  user  model 
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Fig.  4.7.  Uncoded  16k  symbols / sec  QPSK  performance  for  matched  prolate 
spheroidal  windows  with  peak-power-limited  transmitters  and  LOS  users. 


Thus,  PS  windows  of  moderate  length  are  capable  of  handling  both  CT  and 
a particular  form  of  LOS  interference  for  the  closely  packed  FDMA  situation 
outlined  in  Section  3.1.  Since  the  transmitter  and  satellite  windows  need  not 
be  identical,  the  designer  is  free  to  modify  the  transmitter  and  satellite  win- 
dows individually  to  tradeoff  the  different  types  of  interference.  Increasing 
the  length  and  perhaps  reducing  the  cutoff  frequency  of  the  satellite  window 
provides  added  protection  against  stronger  LOS  users  but  at  the  expense  of 
more  ISI.  A longer  transmitter  window  would  help  in  reducing  CT  among  satel- 
lite users  but  also  would  mean  more  of  an  AWGN  energy  loss  under  a peak- power 
limitation  at  the  transmitter  and  more  ISI. 

Hence,  the  close  frequency  packing  of  equiamplitude , unsynchronized,  QPSK 
16  kbps  satellite  users  on  25  kHz  centers  and  in  between  LOS  users  appears  to 

_3 

be  possible  with  only  a small  loss  in  performance  (<  1/2  dB  for  P > 10  ) com- 

s — 

pared  to  the  theoretical  AWGN  case.  With  no  transmitter  power  control  and 

10  dB  interfering  users  the  price  in  performance  could  be  as  large  as  2 or  3 dB 

-3 

above  the  theoretical  for  Pg  2l  10  • As  explained  in  Section  4.5  a few  dB  less 

E,  /N  would  be  possible  with  the  use  of  the  rate  1/2  convolutional  code  men- 
b o 

tioned  previously  to  improve  the  overall  performance  [54,55].  The  simulated 
performance  for  uncoded  BPSK  satellite  users  resulted  in  nearly  the  same  losses 
as  with  QPSK  before  the  coding  improvement. 

Returning  to  the  question  of  maintaining  CE  signaling  by  not  using  special 
transmitter  windows,  simulation  data  for  continuous  phase  FSK  modulation  will 
be  presented  next.  The  emphasis  will  be  on  MSK  rather  than  CPQFSK  since  MSK 
has  a greater  potential  for  controlling  crosstalk  (CT)  among  satellite  users 
at  the  closely  packed  center  frequencies  assumed.  It  was  verified  by  simulation 
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that  CPQFSK  performance  is  indeed  very  poor.  If  a non-rectangular  satellite 
window  is  applied  to  MSK,  it  is  confined  to  two  (binary)  symbol  intervals  to 
avoid  introducing  ISI. 

The  uncoded  bit  error  probability  vs  E^/N q performance  for  10  closely 
packed  MSK  users  at  4800  bps  and  a center  frequency  spacing  of  3750  Hz  is  shown 
in  Fig.  4.8.  In  this  case  a rate  1/2  convolutional  code  is  assumed  so  that 
with  coding  the  information  rate  is  2400  bps.  The  code  is  capable  of  improving 
the  performance  shown  by  approximately  2 or  3 dB  for  the  equiamplitude  user 
case.  Performance  is  adequate  if  there  is  perfect  power  control  of  the  trans- 
mitting users,  but  it  is  evident  that  the  CT  interference  is  too  severe  with- 
out a high  degree  of  power  control.  A satellite  window  here  noticeably  de- 
graded performance. 

If  the  user  spacing  is  tripled  the  precoding  performance  is  very  close  to 
the  theoretical  curve  in  AWGN  even  when  the  interfering  users  are  10  dB  larger 
than  the  user  being  demodulated.  Thus,  the  higher  channel  symbol  rate  MSK  is 
not  very  promising  because  either  a transmitter  power  control  discipline  must 
be  imposed  on  the  users  or  too  much  bandwidth  must  be  sacrificed  which  is  con- 
trary to  the  original  system  objectives.  Conceivably  a rate  3/4  convolutional 
code  and  a 3600  Hz  symbol  rate  would  fare  better  at  the  3750  Hz  spacings  than 
the  1/2  rate  code,  but  it  seems  doubtful  that  coded  MSK  can  beat  uncoded  MSK 
at  these  closely  packed  frequencies. 

Some  very  encouraging  uncoded  MSK  results  are  shown  in  Fig.  4.9.  With 
two  LOS  users  the  performance  is  within  1 dB  of  the  theoretical  bit  error 
probability  curve  even  with  10  dB  interfering  users.  The  application  of  a 
100-point  prolate  spheroidal  satellite  window  has  little  effect  for  interfering 
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Fig.  4.8.  Performance  for  uncoded  4800  bps  minimum  shift  keying  with 
center  frequency  spacing  of  3750  Hz. 
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Fig.  4.9.  Performance  for  uncoded  16k  bps  minimum  shift  keying  with 
LOS  users. 
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users  below  the  10  dB  level.  However,  it  appears  that  interfering  users  as 
large  as  15  dB  above  the  user  being  demodulated  may  be  tolerable  with  a degrada- 
tion of  2 to  5 dB  from  the  ideal  curve  depending  on  whether  or  not  the  window 
is  applied.  Although  the  satellite  window  improves  the  performance  considerably 
in  the  presence  of  20  dB  users,  in  this  case  the  degradation  from  the  ideal 
would  be  4 to  10  dB. 

When  comparing  coded  QPSK  with  uncoded  MSK  on  a bit  error  probability  basis 
one  must  halve  the  ordinates  and  subtract  and  the  coding  gain  (see  Section 

4.5)  from  the  abscissas  of  the  QPSK  curves.  The  results  of  such  a comparison 
are  given  in  Table  1.1. 

Uncoded  MSK  is  attractive  because  peak-power-limited  transmitters  can 
employ  CE  signaling,  there  is  no  ISI  and  hence  no  need  for  compensating  filters 
at  the  satellite,  and  non-rectangular  satellite  windowing  may  be  unnecessary. 

Of  course,  the  convolutional  decoders  employing  the  Viterbi  algorithm  for  the 
rate  1/2  code  would  also  be  absent  from  the  MSK  design.  These  advantages 
should  be  evaluated  carefully  before  selecting  the  MSK  scheme  for  implementa- 
tion, however.  Considering  the  complexity  of  the  demodulation  process  in  either 
QPSK  or  MSK  as  discussed  in  Section  V the  equipment  savings  for  MSK  may  not  be 
very  important.  Furthermore,  coded  QPSK  should  perform  significantly  better 
than  MSK  at  the  smaller  error  probabilities,  e.g.,  P^  = 10  ^ as  indicated  in 
the  next  section. 

4.5  Coding  Gain  for  QPSK 

QPSK  simulations  were  performed  only  for  uncoded  QPSK  but  at  the  coded 
channel  symbol  rate  of  16000  4-ary  symbols  per  second  necessary  for  a 16  kbps 
information  rate  and  a rate  1/2  convolutional  code.  The  actual  coded  QPSK 
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performance  can  be  estimated  from  the  simulated  uncoded  performance  if  the 

theoretical  uncoded  and  coded  performances  are  known.  The  usual  technique 

[55,56]  is  to  take  a coding  degradation  (in  increased  E,/N  required  for  a 

b o 

given  error  probability)  from  the  theoretical  coded  performance  that  is  equal 
to  the  degradation  of  the  simulated  results  from  the  theoretical  uncoded  per- 
formance at  the  same  theoretical  E^/N^  reference  value.  This  procedure  will 
be  illustrated  by  Fig.  4.10  for  some  of  the  QPSK  simulations  of  Fig.  4.7  and 
the  theoretical  performance  of  a rate  1/2  constraint  length  7 convolutional 
code  over  a coherent  AWGN  channel.  For  such  a channel  the  bit  error  probability 
for  convolutionally  coded  QPSK  is  identical  to  that  of  BPSK  [57]  just  as  the 
case  for  uncoded  QPSK  and  BPSK. 

Referring  to  Fig.  4.10  the  theoretical  coded  performance  for  hard  and 

soft  decisions  [42,  54]  indicates  a coding  gain,  i.e.,  the  improvement  in  E./N 

b o 

at  a given  error  probability  over  the  uncoded  performance,  of  approximately  3 
and  5 dB , respectively,  over  the  theoretical  uncoded  performance  for  a bit 
error  probability  of  10  The  actual  coded  performance  for  the  simulation  is 
estimated  as  follows.  For  P^ < 10  ^ and  10  dB  interfering  users,  for  example, 

the  theoretical  soft  decision  coded  E^/N^  is  about  4.8  dB.  At  this  E^/N^  t*le 
degradation  of  the  10  dB  user  simulation  from  the  theoretical  uncoded  curve  is 
approximately  2.0  dB.  Therefore  the  predicted  soft  decision  coding  performance 
for  10  dB  users  is  E^/l^  * 6.8  dB.  By  roughly  extrapolating  the  simulated 
10  dB  user  curve  it  can  be  seen  that  the  actual  coding  gain  could  exceed  the 
theoretical  5 dB.  It  is  generally  true  that  the  actual  coding  gain  for  a 
non-ideal  modem  exceeds  that  predicted  from  the  theoretical  curves.  The  coding 
gain  for  the  hard  decision  case  can  be  estimated  in  a similar  fashion. 
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Fig.  4.10. 


Estimating  coded  QPSK  performance  using  data  of  Fig.  4.7. 
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V.  SATELLITE  IMPLEMENTATION  CONSIDERATIONS 


The  bulk  of  the  satellite  processing  proposed  in  this  report  is  repre- 
sented in  Fig.  5.1.  The  processor  has  four  distinct  components:  1)  a sampler 
or  analog  to  digital  (A/D)  converter;  2)  some  hardware  which  would  window  and 
store  the  data  prior  to  its  demodulation;  3)  a digital  demodulator  which  would 
implement  either  a series  of  discrete  Fourier  transform  (DFT)  calculations  or 
a fast  Fourier  transform  (FFT)  algorithm;  and  4)  a bank  of  estimators  which 
will  be  either  simple  digital  filters  or  Viterbi  decoders.  In  the  first  four 
of  the  six  parts  in  this  section  the  implementation  requirements  for  each  pro- 
cessing component  are  considered  with  respect  to  the  various  system  configurations 
previously  discussed.  In  the  fifth  part  the  number  of  real  multiplications 
required  for  satellite  processing  is  estimated.  In  the  final  part  of  this 
section,  wordlength  requirements  are  considered. 

5.1  The  Sampler 

The  need  for  a sampler  is  common  to  all  of  the  digital  schemes  considered. 

To  process  16  kbps  data  with  50  samples  per  symbol  requires  an  A/D  converter 
clock  rate  of  800  kHz  which  corresponds  to  a sampling  interval  of  1.25  ysec. 
Furthermore  for  both  QPSK  and  MSK  the  received  data  are  complex  and  thus  two 
converters  - one  for  the  real  part  and  one  for  the  imaginary  part  - must  be 
used.  These  schemes  also  put  a slight  extra  burden  on  the  analog  portions  of 
the  system  in  that  the  modulator  and  demodulator  must  each  provide  two  chan- 
nels . 

5.2  The  Windowing  Operation 

Straightforward  MSK  requires  no  windowing  at  all  so  in  that  case  this 
processing  component  may  be  omitted.  However,  satellite  windows  may  be  em- 
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Fig 


5.1 


Satellite  processing  components 


ployed  with  uncoded  MSK  if  strong  interfering  users  are  present.  If  non- 

rectangular  windowing  is  performed  it  must  occur  once  per  symbol  interval, 

and  it  involves  N multiplications.  Since  the  window  is  real  and  the  data 
w 

complex,  this  requires  that  2N  real  multiplications  for  QPSK  (or  MSK)  and 

w 

real  multiplications  for  BPSK  be  performed  in  sec  for  each  user  demodula- 
tion channel.  For  = 100  and  T^=62.5  psec  this  implies  that  a real  multipli- 
cation be  performed  in  312.5  nsec  for  QPSK  (or  MSK)  if  one  multiplier  per  user 
is  allocated. 

Alternately,  20  users  might  be  handled  with  a high  speed  array  multiplier 

capable  of  one  real  multiplication  every  15.625  nsec.  This  may  be  pushing  the 

current  state  of  the  art  [58],  however,  and  would  be  less  reliable  than  using 

several  time-shared  slower  multipliers.  With  synchronized  user  signals  these 

multiplier  requirements  are  eased  by  a factor  U,  the  number  of  users,  since 

N 

w 

then  a single  data  window  can  be  applied  to  all  users  simultaneously.  -j- 
words  of  storage  are  also  required  for  the  window  function  which  is  a real 
and  even  sequence. 

5.3  The  Digital  Demodulator 

Some  form  of  digital  demodulator  is  basic  to  all  of  the  schemes  considered. 
This  component  is  the  most  complicated  part  of  the  system. 

For  QPSK  and  BPSK  a number  of  DFT’s  equal  to  the  number  of  users  must  be 
performed.  If  f^(n+n^)  represents  the  windowed  sequence  for  the  u-th  user 
adjusted  by  the  user  delay  n^,  and  if  F^(o)^)  represents  the  DFT  output  (£-th 
symbol,  u-th  user)  then,  cf.  Eq.  (2.4) 
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where  co^  is  the  actual  radial  frequency  of  the  received  signal  from  the  u-th 
user.  Since  n^  and  co^  are  assumed  known,  Doppler  and  symbol  timing  tracking 
must  be  performed  at  the  receiver. 

The  DFT  of  Eq.  (5.1)  can  be  implemented  using  one  of  the  schemes  shown  in 
Fig.  5.2  for  each  w . In  Fig.  5.2a  the  output  is  sampled  at  the  end  of  each 
symbol  interval.  It  should  also  be  noted  that  these  output  samplers  for  dif- 
ferent  users  will  not  be  synchronized  unless  the  users  are  synchronized  with 
respect  to  symbol  timing. 

The  number  of  real  multiplications  can  be  approximately  halved  if  the 
network  of  Fig.  5.2a  is  replaced  by  the  Goertzel  algorithm  realization  of  Fig. 
5.2b.  The  feedback  operations  are  performed  once  per  input  sample,  but  the 
feedforward  operation  and  the  final  sum  is  performed  only  at  the  end  of  each 
symbol  interval  at  which  time  the  storage  elements  are  zeroed.  For  complex 
data  this  network  must  perform  2N^  + 8 real  multiplications  every  symbol  inter- 
val plus  some  additions.  Interestingly  this  is  roughly  the  same  requirement 
as  for  the  windowing  operation  with  unsynchronized  users.  An  important  dif- 
ference is  that  even  with  synchronized  users  the  discrete  demodulation  must  be 
performed  for  each  user  whereas  the  windowing  only  needs  to  be  performed  once. 

If  the  users  are  synchronized  in  time  and  frequency,  Eq.  (5.1)  can  be 
rewritten  as 
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a.  Straightforward  implementation. 


SAMPLE  TOGETHER 


b.  Realization  requiring  half  the  number  of 
real  multiplications  of  a. 

Fig.  5.2.  Discrete  Fourier  transformer  for  user  u (Goertzel's  algorithm). 
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(5.2) 
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where  fM(n)  is  independent  of  u and  where  represents  the  number  of  equally 
spaced  frequency  channels  allocable.  In  contrast  to  previous  discussion  is 
allowed  to  differ  from  N,  the  number  of  DFT  input  samples.  In  general 
N > (N^  = 64  for  the  simulations)  and  assuming  that  divides  N,  in 

this  case  a faster  implementation  can  be  found  by  rewriting  Eq.  (5.2)  as 
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(5.3) 


This  implementation  requires  a few  more  additions  but  now  requires  only  2N^+4 
real  multiplications  for  complex  data.  The  total  number  of  multiplications  for 
all  users  can  be  approximately  halved  (again,  if  ~ N^/2)  if  users  are  assigned 
paired  frequencies,  i.e.,  if  one  user  has  frequency  k/N^  another  user  should 


N -k 

have  frequency  — - — . If  both  users  are  demodulated  together  the  recursive 
computations  implied  by  Fig.  5.2b  can  be  shared  and  only  the  final  multiplications 
must  be  performed  twice. 

Another  alternative  in  the  synchronized  users  case  is  to  replace  the  bank 
of  individual  channel  demodulators  by  a single  box  which  performs  an  FFT  of 
Nf  points  using  Eq.  (5.3)  and  in  so  doing  demodulates  all  64  channels  (some  of 
which  may  not  be  used)  simultaneously.  It  would  have  to  perform  a 64-point  FFT 
in  62.5  ]Jsec,  which  is  certainly  technically  feasible.  Although  the  number  of 
multiplications  is  not  a totally  reliable  figure  of  merit,  the  FFT  (radix  4) 
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would  require  at  most  384  real  multiplications.  For  QPSK  the  scheme  in  Fig. 

5.2b  modified  for  synchronized,  paired  users  requires  140  real  multiplications 
per  pair  of  users.  The  FFT,  therefore,  becomes  viable  for  6 or  more  synchronized 
users . 

For  16  kbps  users  at  25  kHz  spacings  the  maximum  Doppler  shift  is  about 
±4%  of  the  25  kHz  spacing.  In  this  case  the  implementation  implied  by  Eq.  (5.3) 
may  be  extremely  close  to  that  of  Eq.  (5.1),  which  implies  that  multiplications 
can  be  saved.  Users  cannot  be  paired,  however,  since  this  also  requires  time 
synchronization,  so  additional  savings  cannot  be  realized. 

For  purposes  of  comparison,  an  MSK  demodulator  requires  2N^  real  multi- 
plications per  symbol  interval  per  user  (N  = 50  for  the  uncoded  MSK  simulation) . 

c 

This  computation  could  be  performed  using  an  FFT  if  the  users  were  synchronized 
but  a high  order  FFT  would  be  necessary  and  only  a limited  number  of  its  out- 
puts would  be  used.  The  number  of  multiplications  required  with  an  FFT  implemen- 
tation will  generally  exceed  the  number  required  for  a direct  implementation, 
so  the  FFT  approach  would  probably  not  be  used. 

5.4  The  Estimators 

Each  estimator  reduces  any  intersymbol  interference  (ISI)  and  makes  a 
decision  on  the  received  symbol.  In  the  coded  case  the  Viterbi  algorithm  is 
applied  to  soft  decisions  to  decode  the  transmitted  information.  Since  this 
report  is  not  concerned  with  coding  nothing  further  will  be  said  about  the 
decoder.  Although  ISI  does  not  occur  for  short  (N^  _<  50)  windows  or  for  MSK, 
this  represents  a rather  insignificant  hardware  savings  for  these  two  cases. 

A 2-zero,  1-pole  digital  filter  for  removing  ISI  from  QPSK  requires  eight  real 
multiplications  per  symbol  interval  per  user  plus  some  additions,  an  extremely 


94 


modest  requirement.  If  desired  this  linear  filter  could  be  time  shared  among 
the  users  rather  easily  with  some  buffering  and  control  logic.  For  BPSK  the 
number  of  multiplications  would  be  halved. 

For  the  Viterbi  algorithm  approach  to  combatting  ISI,  4 real  multiplications 
per  symbol  per  user  are  required  for  the  whitening  filter  and  another  32  real 
multiplications  are  required  per  symbol  per  user  by  the  Viterbi  algorithm  it- 
self. Memory  requirements  for  all  of  these  schemes  are  modest. 

5.5  Multiplications  Required 

The  number  of  real  multiplications  required  for  the  satellite  processing 
components  discussed  in  Sections  5.2,  5.3  and  5.4  are  listed  in  Table  5.1. 
Referring  to  this  table,  unwindowed  MSK  takes  1/3  the  multiplications  of  win- 
dowed MSK  and  less  than  1/3  the  multiplications  of  QPSK  for  unsynchronized 
users  with  insignificant  Doppler  shifts.  With  20  synchronized  users  and  insig- 
nificant Doppler  shifts,  QPSK  requires  about  1/3  the  multiplications  of  MSK. 
Unwindowed  MSK  requires  less  than  half  the  multiplications  of  either  windowed 
MSK  or  QPSK  for  unsynchronized  users  when  Doppler  shifts  are  significant. 

Finally,  with  synchronized  users  and  significant  Doppler  shifts,  unwindowed 
MSK  requires  roughly  25%  fewer  multiplications.  Thus,  QPSK  has  fewer  multi- 
plications than  MSK  only  when  Doppler  shifts  are  not  significant  and  when  users 
are  synchronized. 

5.6  Wordlength  Requirements 

Digital  sequences  correspond  to  sequences  of  b-bit  numbers.  Errors  occur 
whenever  numbers  which  require  more  than  b bits  for  their  representation  must 
be  quantized  to  b bits.  Thus  A/D  converters  and  multipliers  have  a certain 
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TABLE  5.1 


NUMBER  OF  REAL  MULTIPLICATIONS  PER  SYMBOL  REQUIRED  FOR  SATELLITE  PROCESSING 

(N  =100,  N =50,  N.  = 64,  U = number  of  satellite  users), 
w c f 


a.  Unsynchronized  Users,  Dopplers  Insignificant 


modulation 

operation 

total 

windowing 

demodulation 

ISI  filter 

QPSK 

200U 

132U 

8U 

340U 

windowed  MSK 

200U 

100U 

0 

300U 

unwindowed  MSK 

0 

100U 

0 

100U 

b.  Synchronized  Users, 

Dopplers  Insignificant 

modulation 

operation 

total 

windowing 

demodulation 

ISI  filter 

QPSK 

200 

384 

8U 

8U  + 584 

windowed  MSK 

200 

100U 

0 

100U  + 200 

unwindowed  MSK 

0 

100U 

0 

100U 

c.  Unsynchronized  Users,  Dopplers  Significant 


modulation 

operation 

total 

windowing 

Dopplers 

demodulation 

ISI  filter 

QPSK 

200U 

208U 

20  8U 

8U 

624U 

windowed 

MSK 

200U 

208U 

100U 

0 

50  8U 

unwindowed 

MSK 

0 

132U 

100U 

0 

232U 

d. 

Synchronized 

Users,  Dopplers 

Significant 

modulation 

operation 

total 

windowing 

Dopplers 

demodulation 

ISI  filter 

QPSK 

200 

208U 

132U 

8U 

348U  + 200 

windowed 

MSK 

200 

208U 

100U 

0 

308U  + 200 

unwindowed 

MSK 

0 

132U 

100U 

0 

232U 
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amount  of  inherent  noise  associated  with  them.  The  energy  represented  by  this 
noise  decreases  exponentially  with  increasing  b.  Thus  by  setting  limits  on  the 
tolerable  noise  a required  wordlength  b can  be  determined.  Hardware  complex- 
ity also  increases  with  b so  the  optimum  b is  the  smallest  one  that  will  intro- 
duce tolerable  error.  Available  hardware  makes  it  particularly  practical  for 
b to  be  a multiple  of  4 bits.  The  noise  analysis  performed  to  obtain  the  follow- 
ing results  is  based  on  a standard  noise  model  [11,59]. 

Assuming  a clip  level  magnitude  of  twice  the  standard  deviation  of 

2 

the  total  A/D  converter  input  of  variance  20^,  the  effect  of  rounding  to 
^A/D  to  ac*d  complex  white  quantization  noise  of  variance 


20 


2 

q 


4(2a2) 

6 


2 2bA/D 


(5. A) 


—2b 

to  the  input,  where  (2  A/D)/12  is  the  variance  of  the  quantization  noise  for 

a real  input  of  unit  magnitude.  The  number  of  quantization  levels  2^A/D 
should  be  chosen  so  that  (5.4)  is  much  less  than  the  variance  of  the  input  with- 
out interfering  users,  namely. 


2 2 2 

2o  = 2o  + A 

s n 


(5.5) 


where  2o^  is  the  variance  of  the  complex  channel  noise  and  A is  the  amplitude 
of  the  user  being  demodulated.  If  U-l  and  E are  the  number  and  relative  power 


level  of  the  interfering  users,  then  the  desired  condition  on  b^^  can  be 
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written  as 


22bA/D  » 


|(l+(l+(U-l)E)^-2 


2a 


JX. 


1 + 


2a 


(5.6) 


2 2 

Replacing  the  user  signal  to  channel  noise  ratio  A /2a  by  (E./N  )(R/B), 

n bo 

with  E./N  = 10,  the  user  data  rate  to  bandwidth  ratio  R/B  = 16/25,  and 
b o 

nine  15  dB  interfering  users,  (5.6)  becomes 
2bA/D  » 12.8 


or 


is  the  required  number  of  bits  in  the  quantization. 


If  white  noise  with  variance  O is  input  to  a DFT  the  output  sequence  is 
2 

white  with  variance  No  , where  N is  the  size  of  the  DFT.  The  effect  of  multi- 
plier quantization  noise  is  equivalent  to  adding  white  noise  at  the  input  of 

1 _2bDFT 

the  DFT  with  variance  — 2 where  b is  the  wordlength  in  the  DFT.  The 

12  Dr  1 

64  ~ZDDFT 

variance  at  the  output  of  a 64-point  DFT  is  therefore  yy  2 . Thus,  if  the 

input  and  output  of  the  DFT  are  to  have  the  same  number  of  significant  fig- 
ures the  internal  wordlength  of  the  DFT  should  contain  an  additional  3 bits. 

One  possible  scheme  then  would  be  to  use  8-bit  A/D  converters.  The  out- 
put of  these  can  then  be  multiplied  by  an  8-bit  window  function  and  input  to 
a 64-point  DFT  with  internal  wordlengths  of  8 bits.  The  output  of  the  DFT 
should  then  be  significant  to  4 bits.  The  wordlength  used  in  the  recursive 
filter  or  the  Viterbi  algorithm  to  remove  ISI  could  be  either  4 bits  or  8 bits. 
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APPENDIX  A 

INDICES  AND  CHARACTERISTICS  OF  SIMULATION  WINDOWS 


WINDOW 

WINDOW 

WINDOW 

WINDOW  SPAN 

CUTOFF 

MAXIMUM  LOSS 

INDEX 

TYPE 

LENGTH 

IN  SYMBOLS 

FREQUENCY 

FOR  PEAK-POWER- 
LIMITED  TRANS- 
MITTER 

N 

w 

NCW 

f (kHz) 
c 

Lt(dB) 

0 

CE 

50 

i 

- 

0 

1 

PS 

64 

3 

12.5 

2.61 

2 

PS 

76 

3 

12.5 

3.07 

3 

PS 

88 

3 

12.5 

3.44 

4 

PS 

100 

3 

12.5 

3.76 

5 

PS 

128 

3 

12.5 

4.36 

6 

DC 

100 

3 

12.5 

3.55 

7 

PS 

50 

1 

12.5 

1.96 

8 

PS 

64 

2 

12.5 

2.61 

9 

PS 

76 

2 

12.5 

3.07 

10 

PS 

88 

2 

12.5 

3.44 

11 

PS 

100 

2 

12.5 

3.76 

12 

DC 

100 

2 

12.5 

3.55 

13 

DC 

150 

3 

12.5 

4.50 

14 

DC 

200 

4 

12.5 

5.22 

15 

PS 

76 

3 

8 

1.89 

16 

PS 

100 

3 

8 

2.61 

17 

PS 

128 

3 

8 

3.26 

18 

PS 

76 

3 

10 

2.47 

19 

PS 

100 

3 

10 

3.20 

20 

PS 

128 

3 

10 

3.83 

21 

DC 

50 

1 

12.5 

15.4 
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APPENDIX  B 

INDICES  AND  CHARACTERISTICS  OF  ISI  FILTERS 


FILTER 

INDEX 

WINDOW  INDEX 
(transmitter/ 
satellite) 

FILTER  SYMBOL 
DELAY 

Nd 

ao 

ai 

a2 

a3 

a4 

E,  /N  +L.  (dB)* 

DOt 

USERS  SPACED 
OR  PACKED 

2 

1/1 

1 

-.1184 

1.685 

-.1190 

- 

- 

.004683 

11.45 

SP 

25 

2/2 

1 

-.2450 

1.902 

-.2496 

- 

- 

.01756 

10 

SP 

26 

3/3 

1 

-.3812 

2.133 

-.3963 

- 

- 

.03671 

10 

SP 

27 

4/4 

1 

-.5200 

2.374 

-.05530 

- 

- 

.05984 

10 

SP 

29 

1/1 

1 

-.1154 

1.664 

-.1163 

- 

- 

.007156 

10 

PK 

30 

2/2 

1 

-.2379 

1.881 

-.2429 

- 

- 

.01889 

10 

PK 

31 

3/3 

1 

-.3801 

2.130 

-.3951 

- 

- 

.03675 

10 

PK 

32 

4/4 

1 

-.5197 

2.373 

-.5527 

- 

- 

.05983 

10 

PK 

35 

5/5 

1 

-.8344 

2.952 

-.9322 

- 

- 

.1104 

10 

PK 

43 

0/13 

1 

-.5240 

2.001 

-.6006 

- 

- 

. 1003 

10 

PK 

46 

2/2 

1 

-.2617 

1.977 

-.2674 

- 

- 

.02071 

13 

PK 

49 

3/3 

1 

-.4231 

2.258 

-.4407 

- 

- 

.04098 

13 

PK 

53 

4/4 

1 

-.5806 

2.531 

-.6195 

- 

- 

.06699 

13 

PK 

54 

4/4 

2 

.1526 

-.6639 

2.565 

-.6700 

.1648 

-.01814 

13 

PK 

58 

5/5 

1 

-.9660 

3.240  - 

-1.087 

- 

- 

.1284 

13 

PK 

64 

0/5 

2 

.0931 

-.4485 

1.860 

-.4514 

.09847 

-.01207 

13 

PK 

84 

1/1 

1 

-.1255 

1.737 

-.1265 

- 

- 

.007689 

13 

PK 

89 

15/15 

1 

-.3286 

1.589 

-.3473 

- 

- 

.05385 

13 

PK 

90 

18/18 

1 

-.3121 

1.775 

-.3244 

- 

- 

.03801 

13 

PK 

93 

12/12 

1 

-.7269 

2.604 

-.7674 

- 

- 

.1005 

13 

PK 

94 

12/12 

2 

.2523 

-.9090 

2.704 

-.8851 

.2593 

-.03463 

13 

PK 

95 

13/13 

1 

-1.315 

3.872  • 

-1.656 

- 

- 

.2061 

13 

PK 

96 

13/13 

2 

.5987 

-2.037 

4.489 

-2.346 

.9005 

-.1190 

13 

PK 

* The  listed  values  are  the  E, 7n  f s for  a rectangular  transmitter  window  (I.  = 0 dB)  so  that  the 

DO  t 

noise  variance  of  each  design  is  independent  of  the  shape  or  length  of  the  transmitter  window. 


MSF 

,085 

,096 

,108 

123 

,096 

,105 

109 

123 

166 

153 

061 

060 

072 

068 

102 

098 

057 

081 

067 

096 

086 

155 

126 


channel  Gaussian 


INPUT 


Fig.  B.l.  MMSE  filter. 
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